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ABSTRACT 


Since 1960 there has been considerable 
interest in experiments involving electron tunnelling 
into superconductors. This new technique, which is 
Simple in principle, provides a direct method of 
observing the modification of the electronic density 
of states that take place in a metal when it becomes 
Superconducting. A forbidden gap forms about the 
Fermi energy in the density of states of the superconductor 
and much smaller variations in the density of states 
occur at energies above the Fermi energy by an amount 
equal to the predominant phonon energies of the metal. 

Investigations of this energy region have been 
made for the metals lead, tin, indium, aluminum, and 
mercury. The results found for lead and tin are in 
agreement with published results, and the results found 
for indium and aluminum are new. Phonon effects should 
be observable in mercury, but specimens are difficult to 
prepare and the results obtained with the one specimen 
tested were uninterpretable. These preliminary results 
are encouraging however, because they hold promise that 
the method used to prepare and mount the mercury specimen 
can yield useable results. 

Values are also reported for the energy gaps 


in the above metals and for the observed amount of 
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broadening of the singularity in the density of states 
at the edge of the energy gap. The possibility of using 


tunnel junctions as thermometers has also been investigated. 
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CHAPTER I 


INTRODUCTION 


A. Motivation 

This thesis reports on a continuation of 
work initiated at the University of Alberta by J. G. 
Adler in 1962. The original motivation for our 
electron tunnelling work has been given by Adler in 
his Doctoral Thesis (1963). 

During the course of Adler's work 
instrumentation was developed which enabled us to 
observe small deviations in the electronic density 
of states of some metals when they became super- 
conducting. The cause of these deviations was not 
entirely clear vat that time, but it. has since »been 
established that they are intimately connected with 
the phonon spectrum of the superconductor. Further 
investigations of this effect are reported here. 

The motivation for the kind of work 
described here has been given by Krumhansel (1962). 
The argument given is that the phonon field is an 
experimentally accessible quantum field with a 
natural high frequency cut-off, and as such it 
Should be an excellent testing ground for quantum 
field theories. 


It should be pointed out that electron- 


a ee 


I ABTIAHO 


to nok%aumtiaes o no adaoget afeelt eld? 
Oo .t os. sired!A lo ysttavevinl eri? ts Seysiviat atow 
tuo tol sobtsvidvom Lenitgiro ei? ,S@0l nt SeLbA 
nf teslbA yd movia aged ead acow antlieanys nortosle 
(€0QL) eised?T Leyovood eid 

j1iow @'sefBA lo seiuoo ef? gulayd 
o¢ eu beldsas doltiw Osqoisves sew notdssromuetere -_ 
¢tLereb otnoudoote oi ni anolttetveb [Lame avisedo - 
-iSdue simtood yedd oodw afestem amos to estage to 
Jon ssw apotjsiveb seed? to sauao sunt ,gatvoubnoo 
ase0 somte as 3f 3ud ,emis tsdz gs teefo yLlerténe 
diviw befoonaos ylstsaatini ers yodt Jedd peratidatee 
seditud .tosoubsootTague add To muadoege nonodq sz 
“ated bestogst sia Josits eint to saotseghtesvat: 
AtoM to baix edj 1% Hoksavtdom sat ‘ia 
.(SO@L) Leensdmutd yd oevia necd ead sted bedker eb 
> 


o 
a 


ns 2i Blatt sotiodg sdt dad’ 2! cevig Ioomugae oat 


ioiw biel? mudmaay¢ sidieassos. ylts ‘sitnem 


a2 done @8 bas ~Tto=duo congo. Ha 


aidisup tol balot gaideet 2i 


~torsaste tans Juo bel 


tunnelling is not the best method of exploring the 

phonon field. With specimens involving polycrystalline 
films, the only type of specimen fabricated to date, one 
can only obtain one characteristic curve, and any 

phonon information must be interpreted from this one 
curve. Peaks in the phonon spectrum cause irregularities 
in the curve, and these irregularities overlap consider- 
ably. Hence the only quantitative information tunnelling 
gives is a few phonon 'critical frequencies',and one does 
not know what symmetry direction in the crystal the 
critical frequencies pertain to. Furthermore we are 
limited by the technique to a few superconductors. This 
is to be compared to neutron spectroscopy, which has 
allowed entire dispersion curves along various symmetry 
directions to be obtained for several metals to a 
precision of about one percent (Brockhouse et. al. 1962, 
Woods et. al. 1962). 

However, neutron spectroscopy has its 
limitations also. The metal studied must be a good 
neutron scatterer, and the neutron source must be intense. 
Even with these two demands satisfied, each experiment 
takes a time of the order of days so that automatic 
recording systems must be used. 

Thus in view of the relative simplicity of 
electron tunnelling work plus the fact that its basic 


mechanism involves an electron phonon interaction at 0°K, 
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it is desirable to get as much information on phonons 
as possible with this method, in spite of the somewhat 
limited nature of the information. 

Even limited data can be of considerable 
value. Consider two phonon calculations based on 
elastic constants and a plausible model, one by Bhatia 
(1955) for sodium and the other by Musgrave (1963)) for 
white tin. The phonon spectrum obtained by Bhatia for 
sodium predicts a specific heat curve which agrees 
with the observed one to within 4 percent. The end point 
omits, spectrum kislrat we. imv.’ The actual end point 
determined by neutron scattering is at 15.8 mv. 
Musgrave's phonon spectrum for tin gives a specific 
heat curve which is 4 percent too high; the calculated 
spectrum is complicated with an end point of 21 mv. The 
actual spectrum as determined by tunnelling is indeed 
cCompudicated with ‘an ‘end point ofl/ #7 mv..-'The vend point 
is unambiguous, and there is no other experimental 
source -orthistinformation at present. 

The calculated phonon spectra mentioned 
above are in error because the authors had only zero 
frequency information available to them (the elastic 
constants) along with integrated information (specific 
heat data). Phillips (1959) has outlined a method for 
calculating phonon spectra if critical point data are 
also available, and these calculations should be far more 


accurate. 
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CHAPTER II 
MICROSCOPIC THEORIES OF SUPERCONDUCTIVITY AS RELATED TO 
ELECTRON TUNNELLING 
A. The BCS Theory 

the [irs successiul microscopic theory of 
superconductivity was evolved by Cooper (1956), and 
Bardeen, Cooper, and Schrieffer (1957 a,b), and is 
usually referred to as the BCS theory. 

While this theory has been superseded in 
elegance by later works, it nonetheless appears to be 
correct in its essential features. It will accordingly 
be taken as the main framework for the descriptions here. 

The BCS theory is based on virtual phonon 
exchange between electron pairs, and is accordingly a 
many body theory. The Hamiltonian Operator + for the 
electron is taken to be the sum of a Bloch term Hy » a 
phonon interaction term H, , and a coulomb repulsion 


term H2 =. Thus in the BCS theory, 


H=HotHit+Ha , (ead) 
where 
Ho= & leal (i-(Ac)) +2 Ce HI 
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in these expressions, Rk is the electron wave 


vector, Rx is the wave vector of an electron at the 
Fermi surface, K is the phonon wave vector, © is the 
Shitiwindex.€ 9 is the Bloch energy measured from the 


Fermi energy, AW is the phonon energy, Myx is the 
matrix element for the electron-phonon interaction 
Calculated for the zero-point amplitude of the lattice 
vibrations, C™(hk,o) creates a Block state aoe a 
C(&,o) annihilates a Bloch state |R,o) » and 


Hike) the single particle number operator is 


x | ; Ik, o> occupied 
Clk) CCR eo) = 
©, |&,o> unoccupied - 
The symbols e@ and have their customary meaning. 
The phonon interaction term is attractive for 
excitation energies len —Egixl < RW, + 
In order that the Gee ce operators in the Hy, term will 
not cause contributions to the energy of alternating sign 
when operating on the various possible wave functions, the 
BCS theory considers only a sub-set of electron states 
cal Raph in which the electrons are paired. The pairing 
ts “such that if “the “state [Re + o> is occupied, the state 
\--R,-Y2> is also. The BCS criterion for super- 
conductivity is then that the attractive interaction term 
H, dominate the Coulomb repulsion term Ha, for this 
sub-set of states. Creation and annihilation operators are 


then defined for the pair states: 
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In order to simplify the mathematics, the 
coefficient of the C ,° C™ operators in Hy, is 
replaced by a negative constant for excitation energies 
in the range |El < “Wp and by zero outside this 
range, where A Wp is a representative virtual phonon 
energy. “ike physical "justification" 6r*this *ts thatthe 
thermal excitation energies in a superconductor will be 
very small compared to the predominant virtual phonon 
energies. The Hy, term of the reduced Hamiltonian 


thus becomes 


Marea: aly g, PR ba 5 lERl<AW, - 


It is worth noting at this stage that because 
GCISthis simplified form for mere » experiments which 
involve excitation energies of the order of AW, 
will probably yield results which disagree with those 
predicted by the BCS theory. This is observed to be the 
case in electron tunnelling experiments, and the regions 
of disagreement are of interest because they are 
connected with the phonon distributions of the metals involved. 
It is not possible to. generalize the BCs theory 
to this wider energy range by simply returning to the 
expression given for Ha however, because lifetime 


effects become important for these highly excited 
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states. The BCS theory has thus been contrived to 
describe the superconducting state for electron 
energies close to the Fermi energy. The results of 
the BCS theory are that 

(a) The single-particle-like excitations have the 


peculiar dispersian law 


En = Ey tub ay 
where ep is the single-particle-like excitation 
energy, ea. is the Bloch energy (both measured from 
the Fermi energy), and A is the real, positive half 
gap energy, which is of course temperature dependent. 
(b) The temperature dependence of A is given by the 


integral equation 


AWp 
oe f dé tanh (@ £224°) 
Seer +} Tanh (24 === ’ 
Nov 2 C24 0% 2 (2.3) 
where WN) is the density of states at the Fermi 


energy for electrons in the normal metal, V is the 
constant appearing in Wi4rea. , and = Ket 5 
where Xp is Boltzmann's constant and ve is the 
temperature. 

With the additional assumption that A 
becomes zero at the transition temperature for the 
Superconductor, it is then possible to tabulate A as 
a function of the reduced temperature ae » which 


Muhlschlegel (1959) has done. 
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An approximate expression for the Muhlischlegel 
. \/ 
tabulation is an =(\- t4)"* . the BCS expression 
Pore) fs 


Ab) = 1-75 Rete . eer 


(c) The electron density of states for the Superconducting 


state is obtained by differentiating equation 2.2: 
pe | 

AdNie) = ANE) GE = Nilo) Rp (Ie 9) ig ) ee 
ole ade cte 
where Re means the ‘real part of'. 

This new density of states is populated 
according to Fermi-Dirac statistics. The electronic 
density of states in a superconductor is thus redistributed 
near the Fermi energy in the manner indicated in Figure 2.1. 

Owing to finite lifetime effects, the singularity 
in the density of states at the edge of the energy gap 
is not found experimentally, but otherwise the BCS theory 


is in excellent accord with experiment for energies within 


a few gap widths of the Fermi energy. 


Bes The Eliashberg Gap Equations and Their Application to 
Tunnelling Results. 


As has been indicated in the preteding section, 
the BCS theory does not take the explicit form of the 
phonon spectrum into account. It is also gauge-variant. 
Bogoliubov et al (1959) have succeeded in reformulating 


the BCS theory in a gauge-invariant manner. Eliashberg 
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Figure 2.1 

The electronic density of states 
near the Fermi energy for a metal 
in the superconducting state, and 
also for the same metal ‘in the 
normal state. The energy origin 
is at the Fermi energy. The small 
energy range does not permit band 


structure to be shown. 
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(1960) has in turn succeeded in extending the 
Bogoliubov theory to include both lifetime effects 
for the excited electrons and an explicit phonon 
distribution. Lifetime effects have been included 
by using the Green's Function (or propagator) 
formalism. In this formalism, energies can be 
complex, so that plane waves decay exponentially 
in time. One thus expects the Eliashberg theory to 
yield an energy gap parameter A= Aw) which is 
a complex function of energy; W is taken as the 
energy variable since the quasi-particle concept of 
the BCS theory is not valid when finite lifetimes are 
taken into account. It is also convenient to measure 
energy from the Fermi energy. 

Throughout the remainder of this thesis, 
Planck's constant #A will be taken equal to unity, 
as will the electron charge €@ and Boltzmann's 
constant 2 - this is\Gommon practice in theoretical 
Dapers, and 1t causes no difficulty in the experimental 
work to be described here in that the connection between 
natural frequencies, electron energies due to applied 
voltages, and thermal energies is readily obtained by 
reinserting the appropriate constants. That is, W=eVeT 
becomes AW=ev=RaT for purposes of calculation. 
Kittel (1960) has given an excellent tabulation of the 


various ratios of the physical constants which facilitates 
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such conversions. The ratio of most interest here is 
2/24 = |1-60 °K/mv 

Returning to the Eliashberg theory, it is indeed 
found that A is a complex function of WW) =, but the 
equations which determine A are not very tractable. 
In the Eliashberg theory ,A(W) is given in terms of 
a gap function PD (v0) and a renormalization function 

Z(W) according to the equation 


Ale) = PMV/z (uy 


The equations which determine Q and £ are 


Wee =\/p 
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where gw) is the normalized phonon spectrum, o(W) 

is an effective electron-phonon coupling parameter, U 

is a Coulomb pseudo-potential defined to include only 

the interactions outside a band of energies |WI| < Wee 

which is large compared to the Debye energy, Ao= Paes in.y 


is the gap parameter at the edge of the energy gap, 
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and Duy 0) = ( K+ Gor bot) is one of the phonon 
propagator terms. ) atid? «or are assumed to be 
isotropic and homogeneous. 

Once the function A(W) has been determined, 
it only remains for-our purposes toefind\out»how it 
affects the single-particle-like density of states ina 
Superconductor. Schrieffer et al.(1963) have been able 
to show, again with the aid of the propagator formalisn, 


that the effective tunnelling density of states ina 


Superconductor, AN = NW) 3 tSeeives by 
AW 
ae eg) 
= yl EE), ‘ 
Nd) = NiO Re | | (Acca | | ? 


where as before N60) is the density of Bloch states at 
the Fermi energy for a normal metal. For energies where 
phonon effects are important, W>yla\ w equation (2a*7 ) 
may be expanded to first order in ee 
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where A(W) = A ,(W) + ( LA 2() ) 
and “os (—\)e2 . 


It is perhaps desirable at this stage to show 


briefly how electron tunnelling results are connected with 


NrW) | 


the normalized density of states for a superconductor, Wal 
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metallic films and a dielectric layer arranged in the 
geometry of a parallel plate capacitor. The dielectric 
layer is sufficiently thin that electrons may pass from 
one metallic film to the other by 'tunnelling' through 
the potential barrier provided by the dielectric layer. 
The tunnelling process is that predicted by quantum 
mechanics, and it is a temperature independent process. 

A simple model for the tunnel current exists 
for the case of an n=-s system* at absolute zero which 
exhibits negligible band structure and also has a 
constant barrier transition probability P for energies 
near the Fermi energy We 5 

This model is shown in Figure 2.2 along with 
the current-voltage (i-v) diagram and its derivative. 

The subscripts a and b refer to the two metals involved in 
the tunnel junction. 

The current-voltage (i-v) expression for this 


case is seen immediately from the diagram to be 


4 
LN) = Const: p/ N-(W)o NO), AW 
(ey 
The derivative of this expression is 


ou 


| = const» P- NW NO) 
Av Is 


* i.e., a tunnel junction in which one metal is in the 
normal state and the other is in the superconducting state. 


an3 nt begnaetis ieysl str stootetDve' bas emt ft obte 
as in 
ntttosrfelb SAT .toittosqso esaiqy feflsisg 6 ‘to pes 
mort easqg yam enorroels ftaiz aba yisastorrive eho 
auoidd 'gnkifennyt' yd istide sig cpl eacaanak 
~teyal obstoefstd os yo bebtvotq <sittad ftettnastoq en? 
@utasup yd bavotbe rq shad Bt Sseponyg gutifonnaus 6dT 


gee001g fnebasqebal erusstsqmed & gt 22° be .eoinedo 


¥ 
= « 
s - 


stetxe tne rio fennud ans nol fosom siqute A —— 
dotdw oven saduloeds Js imataye 2=f1 kh to en 
est oeis bas siutouite baad sidigthigen sF 
asitatens ‘1o'}! q yoiiltdsdotg nolyLerens ies on 
; 70) yuxene Lost edt 
ddtw anols &.S saugtt at pwode ef Pepom efar * 
_ ovitdsviveb 231 baa msigsi6 (v-2) egatlov-s3noTiwe 
at bsvlovnl efesem ows sdz oy a9'lsy d Das atqtroadue sa? 
Liiokdoaut, Lennud etd 

7 


- etddg i912 nobeesiqus. (v—t) susilov-dassaus eat . ier * 


ae 


6¢ o3 matresib edd mort i head nes @ 


ee -tenod = vi 


st nolesstgxs onde to oF td 


HOW +l) «Fad 


a , = “ 


oAs af at Lsisom iro. Bokitw rt 
etete gatscvbnoorsque ong nk at, a9 


14. 


Figure 2.2 


A model for the tunnel current 

ang tne électrical ¢haracteristics 
resulting from this model. The 
distribution Nyr(w-v)p is 
given by 
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fe I 
The subscripts has been appended to indicate that one of 
the metals involved is in the superconducting state. 


Repeating this analysis for the case where both metals 


are normal yields 
at 
Av !n 


The ratio of these two derivatives, o-lv) , is thus 


See narT = Te re Nola” 


ow) = de /s) = Nylvjo 

Vis / AVIN Nlole 
so that in this approximation the normalized dynamic 
conductance —“¢ fora’ Tunnels juncté on- ds equivalent 
to the normalized density of states in the superconductor. 
A more detailed analysis of this equivalence will be 
given in later sections. 

Having thus shown that Niki) is’ a directly 
measurable quantity, we now return to equation 2.8. It is 
immediately clear that the equations 2.6 must be solved if 
equation 2.8 is to be compared with experiment. 

Schrieffer et al. have solved the Eliashberg 
equations 2.6 numerically for a model chosen to represent 
lead. In this model the parameters X(w) and U- were 
taken to be constants and 3) was taken to be the sum 
of two Lorentzian-shaped peaks* centered at energies 


corresponding to the predominant transverse and longitudinal 


* A derivation of the Lorentz line-width expression is 
given in Dicke and Wittke (1960). 
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phonon energies. The phonon model distributions are 
shown in Figure 5.3 along with our observed OO results 
for lead. it may be seen from Figure 5.3 that the & 
characteristic shows a pronounced decrease at the phonon 
energies measured from the gap edge singularity. This 
experimental behaviour is in excellent agreement with 
Che machine! calculations, of Schrieffer.et.al., Their 
computations show a sharp increase in Lv.) and a less 
rapid decrease in A&W) at an energy just beyond each 
peak in the phonon distribution. (i.e., at an energy 

V= ey i often, ) where PM is slightly larger 
than the phonon peak energy). 

At these energies it will be noticed that 
excited electron states may readily decay by phonon 
emission to the states available at the gap edge sing- 
ularity. Thus as the voltage applied across a tunnel 
junction is increased past a predominant phonon energy, 
the decay rate for the excited electrons being injected 


intovuthe metal increases ,and) this in«turn, causes a 


decrease in the © characteristic due to the increase 
of the,.Ag term in the density of states 
expression. 


An alternative method of solving the Eliashberg 
equations 2.6 has been given by Scalapino and Anderson (1964). 
These authors have obtained approximate solutions of 


equations 2.6 by using a re-iterative procedure based on 
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Singularities. The first iteration consisted of using 
an Einstein phonon distribution G CW Wo) = 6 (W)=),) 
in the gap equations 2.6. The resulting functions 
QD ( W, Wo) and ZRCwNGsa4) were found to 
be singular at W=W>5 =. The second iteration then 
consisted of using the singular parts of the Einstein 
results Q(W,We) » £ (Wd, We) in equations 2.6 
along with a more plausible phonon distribution in order 
to calculate new values of (Q and # . The forms 
of phonon distribution used were those predicted by 
Van Hove (1953) for phonon energies close to the critical 
point energies Ke - The critical point energies are 
those at which the (0 (q) surfaces for the phonons 
have a maximum, a minimum, or a saddle point. (The 
d)(G) Surfaces are just the three dimensional version 
of the usual phonon dispersion curve K)(q) versus GQ for 
a given direction of the phonon wave vector q a3 
Scalapino and Anderson found that certain of 
the Van Hove forms for Qiu) gave yY and # #£=functions 
in their second iteration that produced logarithmic 
Singularities in the derivative AS of the @ character- 
istic*, at voltages V=WetAo - These authors state 


that their second iteration results are exact as to the form 


om 


=. kenge the: C characteristic for the model given earlier 
in this section. 
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of the singularities and nearly so as to their magnitude. 
One of the Scalapino and Anderson results of 
particular interest to us is that at a voltage corresponding 


to the end point of the phonon spectrum one should find 


limit Aone = BC In |Wee nevi 

V-> Ag+tic dv Wt { 219) 
and 

limit ( ds + 3S: ae TGc 

oo div DotWete AV 1A G+ O.-€ UO ¢ 
where C is a constant, 3 and ¥ are slowly varying 
functions related to O)¢ Pee ee rs is the predominant 
transverse phonon energy, AX05 is the energy gap at 


absolute zero, measured at the edge of the gap for the 
Superconducting member of the n=s tunnel junction, and 

We is the end point energy for the phonon spectrum 
of the superconductor involved. 

Equations 2.9 may be viewed as a quantitative 
Gescription of the behaviour of the experimentally obtained 
GC: characteristics for applied voltages corresponding 

to the end point of the phonon spectrum. If there is a 
logarithmic singularity in the He characteristic near 
the end point of the phonon spectrum, then one can determine 
the end point energy unambiguously. It is unfortunate that 
the i characteristics currently being obtained 


experimentally do not exhibit well-developed singularities, 
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so that the Scalapino and Anderson analysis may not be ex- 
ploited to the extent that one might desire. 

It is of course impossible for the observed 
characteristics to be truly singular when the observed ao 
characteristics do not exhibit well developed singularities 
at the edge of the energy gap. 

It is because of this lack of well developed sing- 
ularities that we limit our analysis to the end point of the 
total phonon spectrum if further data is not available. One 
might, for example, expect singularities of the form (2.9) 
to occur at the end point of each polarization branch of the 
phonon spectrum. Such a result is not guaranteed however 
because the effect of a saddle point in one branch of the phonon 
dispersion surface may cancell the effect of a maximum in 
another branch (Van Hove). Conversely, the observation of a 
singularity of the form (2.9) in observed as character- 
istics at an energy corresponding to an intermediate point in 
the total phonon spectrum need not indicate the end point of 
one of the polarization branches, since such a result could 
equally well be due to a saddle point. The problem is further 
complicated by the possibility of surface maxima which do not 
admit a quadratic term in their Taylor expansions. 

The only completely general result which we have for our 


tunnelling results then is equation (2.9) for the end point of 


the total phonon spectrum. 
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If information is available from other sources 
however one may be able to make more extensive critical point 
analysis of tunnelling data. An example of this possibility 
may be seen in the results obtained for aluminum (Figure 5.11). 
In this case one has Walker's x-ray results for reference, and 
the end points of the three polarization branches may be 
identified in the tunnelling results. The tunnelling results 
thus complement the x-ray results very nicely, in that the 
energy scale for the x-ray results is the most uncertain 
feature of that method, whereas it is the most certain feature 


of the tunnelling method. 
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CHAPTER III 


ANALYSIS OF ELECTRON TUNNELLING RESULTS 


In the preceding section a tunnel integral was 
Solved for a restricted case. We now wish to obtain 
Solutions to tunnel integrals which are subject to fewer 
FestraculLons:, 

Harrison (1961) has shown that band structure 
effects should not be observed in electron tunnelling 
characteristics for two normal metals. The argument 
given by Harrison is based on an independent particle 
model, and he finds that density-of-states terms 
(proportional to Vg » where V9 is the group velocity) 
cancel current terms (proportional to Va | ) in the 
tunnel current integral. The density of states 
distributions which we observe with electron tunnelling 
into superconductors arise from the subtleties of the 
many-particle system, According to the arguments of 
Bardeen (1961), and Harrison, the density of states 
observed with electron tunnelling into superconductors 
is that of the quasi-particles only. When calculating 
tunnel current integrals, we may use models that are 
free of band structure, and replace the density of states 
in the normal metal, N(w)aq , by the density of states 
at the Fermi surface, WN (0) ; Obviously, interesting 


results will only be obtained when at least one of the 


metals is superconducting. 
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Two such cases are illustrated in Figures 3.1 and 
3.2. The first line of these figures depicts the 
normalized tunnelling density of states = Nrieils as 
functions of energy. The subscripts §$ and n- refer to 
the superconducting and normal states respectively, and 
A= 

The distributions Ps are populated according to 
Fermi-Dirac statistics. The number of thermally excited 
electrons in aluminum is negligibly small at the lowest 
temperatures available in the present apparatus, so, by 
varying the temperature, the number of thermally 
excited electrons in the aluminum may be changed at will. 

Figure 3.1 shows some thermally excited electrons, 
and it may be seen that the small peak which they produce 
in the current-voltage (i-v) characteristic provides a 
means of determining the values of the energy gap parameters 

Da. and Ay, ° 

It is possible in certain cases to relate the 
distribution Pes. to the observed i-v characteristics 
without knowing the explicit distribution Psy s rer 
example, consider the situation shown in Figure 3.e. Here 
there are no thermally excited electrons and the Ps, 
distribution is taken to be free of phonon structure. Both 
these conditions may be met quite well with a sufficiently 
low temperature and a suitable choice of metals. The Fermi 


function Akin may accordingly be replaced by a unit 
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Figure 3.1 
ITllustration of the use of 
thermally excited electrons 
for determining energy gap 


values. 
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Figure 3.¢e 
Illustration of a method of observing 
phonon structure effects in tunnelling 


characteristics. 
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step function: 


The expression for the current is readily seen 


from the diagram* to be 


c | Plwv) Row) Ry (w-v) {8 (w-v) ~ O(w) }dur, 


where P(w,V) is the probability of an electron crossing 
the potential barrier, and C is a constant. Upon taking 
the derivative of this expression with respect to voltage, 


and noting that RP.) =O , one obtains 


+S SS = f 2B) R&) Ri (o-v) dos 


Vv 
+ [ Plow) R bw) 0RelW-v) dw 
6 


OV 
I,e+la ro 
Noting further that i) IR CW-Y) dw = 4 and 
° OV 
that 0p (w-v) is@sero*’except; near #G) My ANG » 


OV 


one obtains 
pe P(v-Ab,v) Ro v-Ab) 


If a similar analysis is performed when two normal 


* or, alternatively, see Shapiro et al (1962). 
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metals are involved it yields 
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C davinn f oN 
de to AL 
Denoting the ratio of aa |e a 
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an Vv 
OP er 
P(v,v) + i ae AW (3-1) 
Since AAD eax may be quite small (of the order of 
1032) snd stall. be of anterest, ibeis nora all glear that 


P (wW,V) may be treated as a constant in this 


expression. Information about the behaviour of P(W,v) may 
dé 

av inn 
P(W,Vv) about the origin yields 


be obtained from For example, expanding 
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Figure 3.3 shows this ratio for an A1-~A1,03-In tunnel 
Junction and it may be seen that the neglect of 

higher order terms in the above expression is not serious 
Since the curve is reasonably linear. In this particular 


case one obtains 


pakes s oP =- 2x (agro 
Plo,o) OV 10,0 
Since Figure 3.3 is reasonably linear, one may 
replace Urs WV) by the constant oP in the 
OV 0,0 


integrals in equation 3.1. The integral in the numerator 
then becomes negligible* as does the integral in the 
denominator. Expanding P(v,v-a,) about P(v,v) then 


yields 
O-| = RP WHAy -i 


For the case being considered (Al-In), Ag\=0-2mv , 


so that 
F =5 
Sete ee von) = 1 = 410 


The contribution due to the variation of the 
potential barrier is thus negligible since Oo may 
only be measured to a few parts in ten thousand with 
the present apparatus. 


The major limitation on the expression 


* The energy states which in the normal metal would be in the 
energy gap region form the gap edge singularity in the 
distributions, and M=1 Forirtw 7A since phonon 
effects are generally small. 


pen 


; ae 


fennut al«pUgih=la as OT Obdet, aidd awoiule sD 
DP A 


eT 
Yo tosinen ody tadd sese od yam 32 Ons MOLIOND] 
es 
15 ‘5 


uotues Jon at agheéetqxs avods ens nt anred tabtte , aa 


¢ 


. 
5 a 
@ 


ssfwottasg etdt nl .teeakti yldsnoeset 3 onsite ite 


inte 


- aitstdo ono 28) aia 


vent smo ,1sentl eldsnoeses af £.f omplt eaaks 


. 
aK | | ; 
‘edt at | oe B tdsgeroo saz Yd (vy 6 b} 


uz — 


: 6.0) VO v6 


- aa an 
1nofdetemun oft al Isagsdnt eat .i.f vo tdaupa ab ehetgesat 
ae ¢ 7 ” 


: - 2 a an 
. sd7 at {suednl oft seod as *sldtgtiges asmooed node 
0 : rar on -_ 
ian (yw Ipods (ygd-V Vv? antbnas@x3S ,tovanimoned 
| G6 7 
| so = — }— (sa=-¥v) - 2 


7 
vin S-O =), .(ai-1A) benepisnee gnied sano sad 10% 


isdv wes “8 


“Oink = b- (yaaw) a+ Ps 


sit to noltgstusv Savy 63 Sip oO Lsudl epg se 2 eo 
x i we, 
lo ie i 


vem SD sonte eidteiigen aude at aoiliid i fat 


43itw basevodt ast nt edasq’ wet s ot bowie +: 


oe ie ms ‘nd (asad, is A 
| ot Ves 
ra aoa oes 


. 


q 
=e ? 
: ? 


Figure 4.3 
Transition probability test for 


specimen In-38 at 4.6°K 
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eh, = R,(v- Ae) for the situation depicted in 
Figure 3.2 is in the assumption about ee 
being localized. It is only valid if this derivative is 
non-zero over a region small compared to that occupied 
by the variations being observed in Ro . 

Distributions performing the function of Pe 

in the foregoing example will be called probe distributions. 
It is readily apparent that either a superconducting or a 
normal probe distribution at non-zero temperatures will 
Gistort the observed distribution. The superconducting 
probe distribution tends to exaggerate variations in the 
observed distribution while the normal probe distribution 
tends to diminish them. For this reason superconducting 
probe distributions have frequently been used for this work. 


An attempt to evaluate the distortion which they introduce 


is given in Appendix IV. 
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CHAPTER IV 


EXPERIMENTAL METHOD 


A. Specimen Preparation and Mounting 


Tunnel junction specimens have been prepared, and 
mounted in the He3 cryostat, in the manner described in detail 
by Adler (1963). This method of specimen preparation results 
in a four terminal tunnel junction on a glass substrate. The 
prepared films have been from 300 to 600 A thick, and have been 
polycrystalline, the crystallites being less than 3000 A in 
size according to x-ray transmission photographs.* The tin 
films have shown some preferred orientation of crystallites. 
The substrate is enclosed in a small, gas-tight specimen 
chamber containing one atmosphere of helium gas at room 
temperature, and this chamber is attached to the He? chamber 
of the cryostat with a 0.040" diameter copper wire. Electrical 
leads to the specimen chamber are of lead coated 0.0025" 
manganim wire. A 220 ohm Speer carbon resistor is mounted on 
the He3 chamber to serve as a thermometer. 

Some particular aspects of specimen preparation will 
now be discussed. It will be convenient to designate the first 
metallic film deposited on the glass substrate, the base layer; 
the dielectric layer formed on the base layer, the barrier 
layer; and the metallic film deposited across the barriean 
layer, the cover layer. The cover layer will sometimes be 
designated Ma , when it may be any metal. 
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* The author is indebted to Dr. S. S. Sheinin for the x-ray results 
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musvihe, base Layer) Metal 
Al-A1303-Ma systems have been the only type 


studied here. The reasons for choosing aluminum as the 
Base “layer“are that a-film,of, it is .easily..evaporated, its 
Surface can readily be oxidized to form the barrier layer, 
and it has been studied extensively by other workers. 
Aluminum also becomes superconducting at 1.20°K, so that 
aluminum base layers give a superconducting probe 
distribution at lower temperatures. This is an advantage 
for preliminary phonon work in that such a distribution 
tends to exaggerate any phonon effects that may be present. 
It 18 a disadvantage for critical point analysis, but perhaps 
not a serious one. 

Of the other metals available for a base layer, 
tin appears, from the literature, to be the next most suitable. 
However, the transition temperature of tin is 3.73°K, so 
that it offers an advantage over aluminum as a base layer 
for the observation of phonon effects only if it is also 
the metal being examined. 

Magnesium has been used as a normal base layer by 
Giaever. It is not as readily handled in an evaporator 
because the magnesium atoms do not appear to stick readily 
to anything but the very cleanest surfaces. Greater care 
must be exercised in masking and surface preparation 
procedures than with aluminum, Also magnesium oxide in 


bulk has about one half the volume of the magnesium from 
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which it was formed (Holland (1961) , page 45\ ), so that 
magnesium oxide layers might be porous and contain 
filamentary bridges connecting the cover metal to the base. 
Such bridges of the superconducting metal may completely mask 
the tunnel current but weak filamentary bridges of 

magnesium might be quite tolerable since they would 

remain normal at all temperatures. Magnesium could perhaps 
prove to be a very useful base layer for the observation of 


phonon effects. 


CVS CiueS barrier Layer 
An extensive study of Al-A1503-Al tunnel junctions 


has been made by Fisher and Giaever (1961). The most 
striking feature of their results is the large scatter on 
their graphs of tunnel resistance as a function of barrier 
thickness. The barrier thicknesses were determined by 
capacitance measurements. The observed scatter of about 
20% in the thickness (or the capacitance) could be 

readily explained by an oxide layer of non-uniform thickness, 
Since the capacitance depends on the average of the 
reciprocal of the barrier thickness while the tunnel 
resistance depends on the average of an exponential of 

the barrier thickness. A non-uniform oxide layer would 
also be consistent with the filamentary bridges which are 
often observed. 


Capacitance measurements would thus appear to be 
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a poor way of determining the effective tunnel thickness 

of a barrier layer, but they do at least give some 
additional information when used in conjunction with 

the tunnel resistance measurements. A small R-C bridge 

has accordingly been devised which is suitable for 
measuring both the resistance and the capacitance of 

our tunnel junctions. This bridge is described in 

Appendix II. It covers a larger dissipation range than 
commercial bridges and cannot subject the tunnel junction 
to more than 50 mv. The first feature is desirable 

because of the strong exponential dependence of the 

tunnel resistance on the barrier thickness, while the 
second is necessary because voltages of the order of one 
volt will exceed the dielectric strength of the thin 
barrier layer and destroy the tunnel junction. In addition 
to these features, a display mode is used which allows 
small deviations from a linear current-voltage characteristic 
to be observed. The major limitation of the bridge circuit 
is that it only provides a two terminal measurement, so that 
for accurate results the film resistance must be negligible 
compared to the tunnel resistance, 

It is possible to make four terminal resistance 
measurements using the curve tracer described in Appendix I. 
The display available with this apparatus can be quite 
sensitive, so that the voltage dependence of the tunnel 


resistance may be clearly exhibited. 
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Unfortunately, specimens with weak filamentary 
bridges will appear good when subjected to any of the above 
measurements at temperatures where both metals are normal 
conductors. The final test of a prepared specimen is thus 
only obtained after it has been mounted in the cryostat and 
cooled down to about 0.5°K. At this temperature any 
filamentary bridges of aluminum will be superconducting 
and their effect will be readily apparent at the origin 


of the current-voltage characteristic. 


Den parrier Layer Formation 


A trial-and-error procedure has been adopted 
for oxidizing the surface of the aluminum base layer 
because: 
1. the effective barrier thickness is strongly 
dependent on the metal of the cover layer 
(Handy 1962), 
2. .the thickness. of.the .oxide layer.on ,.the 
aluminum appears to be dependent upon 
the moisture content of the air, and 
3. the oxidation time and temperature do not 
appear to be critical. 
The third reason may appear surprising, particularly when 
these are the variables used by experimentalists to control 
the barrier layers they obtain. Our results indicate that 


an Alj03 layer approaches a saturation thickness that is 
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very near that required for good tunnel junctions. Controlling 
the time and temperature to within ten percent seems to 

give reproducible results for a given cover layer on a 

Given day. The results vary from day to day in a way that 
appears to be associated with the humidity of the Saris: 

however attempts to use deliberately moistened air 

produced more consistent barrier layers but ones through 

which temperature dependent (i.e., non-tunnelling) 

conduction processes were large. Further work might 

develop a consistent procedure that would avoid the 


filamentary bridges that we sometimes obtain. 


BE. Alternative Barrier Layers _ 


Shapiro et al (1962) have reported successful 
tunnel junctions using barium stearate as the dielectric 
layer. They have since indicated* that such barriers have 
a tendency to rupture at low temperatures. 

We have tried evaporating SiO and MgFo dielectric 
layers but have not been able to obtain tunnel junctions in 
this manner. It is quite easy to prepare capacitors in this 
way, but as the dielectric layer is decreased in thickness 
below a few hundred A (as given by the capacitance), short 


newt , 
circuits appear. Since layers about 50 A thick are 


* private communication 
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required before appreciable electron tunnelling can occur 
it has not been possible to use these materials as barrier 


layers. 


F. Lead Connections 

The method, described by Adler, of soldering 
leads to the substrate before the tunnel junction is made 
requires an oxidation temperature lower than 140°C so that 
the indium solder will not melt. The surface of the indium 
contact must be scraped off before the film is deposited 
in order to get an oxide-free indium interface. 

The connections were made somewhat differently for 
specimen Sn-41,. Copper contacts were deposited on the 
Substrate by evaporation ahead of the base layer. Leads 
were then indium soldered to the copper areas after the 
specimen was completed. This procedure allowed us to use 
a higher oxidation temperature (220°C). The specimen which 
resulted was free of filamentary conduction at 0.5°K, but 
this could be a fortuitous result. 

Several attempts were made to prepare mercury 
specimens. A different arrangement was adopted so that 
the mercury could be deposited in the evaporator onto a cold 
barrier layer ( T= 80°K ) and be kept cold thereafter. It 
was found that unusually thick oxide layers were necessary 
in order to get a useable specimen. This is rather 
surprising in view of Handy's (1962) work in which the cover 


layer atoms were found to penetrate the oxide layer in 
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accordance with their atomic diameters. 
Only one specimen was obtained in this way which 
appeared useable at 80°K, but when it was mounted in the 


cryostat, it showed excessive filamentary conduction at 4°K. 


G. Thermometer Calibration 

Temperatures quoted in this thesis are those of 
the He3 chamber inferred from the resistance of a 220 ohm 
'Speer' carbon resistor. The experimental arrangement 
has been described by Adler (1963). Adler's method of 
calibrating the resistance thermometer consisted of using 
the He3 gas as a thermal link between the He3 and He’ 
chambers of the cryostat. The vapour pressure of the 
liquid in the He4 chamber then gave the temperature of 
the resistor. This method of calibration requires that 
the He? chamber be isothermal with the He4 chamber. 

As an alternative to this method we mounted pure 
aluminum and cadmium wires on the He? chamber in order to 
get calibration values for the resistor at 1.20°5. and 
0.52°K where the resistance of the wires dissappeared. 
Another calibration value was obtained at the boiling 
temperature of the helium in the dewar with exchange gas 
in the vacuum jacket of the cryostat. The calibration 
points obtained in this way should be reasonably accurate, 


yet using these points in the Clement-Quinnell (1952) 
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produced resistance-temperature relations similar to the 
one shown in Figure 4.1. 

ihe reversal of the interpolation curve at 7.3°K 
Suggests that the interpolation formula is not well suited 
to this resistor over this temperature range. A further 
calibration value at 80°K has enabled a visual estimate 
of a more plausible resistance-temperature behaviour to 
be made (see Figure 4.1). From this estimate we assign an 
uncertainty of 0.2°K to our interpolated temperatures in 
the 4°K to 1°K range, the uncertainty diminishing near 
calibration points. It will be seen in the following 
chapter that this rather large uncertainty in temperature 
does not affect our results due to the nature of the 
experiments. 

It has also been possible to directly determine 
the temperature of Al-A1503-Al tunnel junctions at low 
temperatures (0.6 to 0.3°K). The current observed at 
applied voltages V< 24 for such a tunnel junction is 
due to electrons which are thermally excited across the 
energy gap. In principle at least, its magnitude gives 
an absolute measurement of the temperature of the junction 
since the width of the energy gap is also in direct 
evidence on the current-voltage characteristics. Preliminary 
calculations of this nature are discussed in Appendix III. 
The results obtained are in agreement with the thermometer 


to within 0.1°K in the 0.3 to 0.6°K temperature range. 
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Figure 4.1 


Interpolation results for a 220 


carbon»resistance thermometer. 
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As before, this uncertainty in temperature does not affect 
our experimental results. The result is interesting 
however because it indicates that the active part of the 
tunnel junction is slightly warmer than the He3 chamber. 
The author believes that the temperature gradient occurs 
mainly between the specimen chamber and the He chamber, 
indicating that 
a) the vacuum chamber has not been evacuated 
sufficiently, and 
b) that the specimen chamber has not been 
thermally anchored to the He3 chamber as 
firmly as one might desire. 
More reliable methods of temperature measurement 
Should shortly be available in this laboratory through the 
work of Dr. F. D. Manchester and Miss J. Gould, who are 
currently investigating the properties of semiconducting 
resistance thermometers. A major source of difficulty with 
carbon resistance thermometers is that their characteristics 
change with thermal cycling between 4°K and room temperature. 


Semiconducting thermometers are not expected to do this. 


H. Instrumentation 

The phonon structure present in © characteristics 
is usually very weak, so that it is desirable to observe the 
dynamic conductance oe in as sensitive a manner as 
possible. Appendix I is a reprint of an article which 


describes apparatus developed for this purpose. Since the 
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di 


dt characteristics are time independent, any assessment 
of equipment sensitivity is somewhat arbitrary, It appears 
however that the noise level of the detecting amplifier 
affords a reasonable basis for comparing systems which 
utilize bridge circuits, and on this basis the system 
described here is competitive with the system described 
by Thomas and Klein (1963) and used by Rowell et al 
(1963) in similar work. 
It is worth noting that the instrumentation is 

not the only factor which determines the uncertainty in 
the observed © characteristics. Any non-tunnelling 
conduction present in the specimen will usually be erratic, 
and this will also cause an uncertainty in the observed 

CG values. It is thus possible to have a ‘noisy’ 
specimen, and a rather ideal specimen is required if the 
specimen noise is to be much smaller than the noise 


level of our present instrumentation. 
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CHAPTER V 


EXPERIMENTAL PHONON STRUCTURE RESULTS 


A. Probe Distribution Details 

The base layer metal used for all specimens 
Studied to date has been aluminum, so that we wish to 
justify the use of the density of states distribution 
of aluminum as a probe distribution. 

The, observed, © characteristic for an Al-Al 
tunnel junction is shown in Figure 5.11. The coordinates 
for this figure have been chosen to emphasize phonon 
effects. The effects in evidence are twice as strong 
as they would be for an n-s tunnel junction, so that 
it may be seen that the use of a superconducting aluminum 
probe distribution has contributed a negligible amount of 
phonon structure to the O characteristics observed 
for lead, tin, indium, and mercury. These metals also 
have larger energy gap values than aluminum, so that the 
density of states distribution of aluminum in the 
superconducting state satisfies our probe distribution 
requirements for these metals reasonably well. 

It will sometimes be convenient to refer to the 
normalized dynamic conductance characteristic, 

oj 30]. / Sl, » as the first derivative of the 
i-v characteristic of the tunnel junction when one or both 


of the metals involved is superconducting. Similarly, the 
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derivative of the G characteristic, AZ will 

AV 9 
Sometimes be referred to as the second derivative 
characteristic. 

In order that phonon structure energies may be 
read directly from the energy coordinate, the figures 
emphasizing phonon structure will have their energy 
origin at Apt Ama » where Ma is the cover 
layer metal involved. 

Figure 5.1 shows the current-voltage (i-v) 
characteristic and the first derivative characteristic 
in the energy gap region for specimen Al-42. The energy 
half-gap value has been determined by extrapolating 
the linear part of the 'current jump' on the i-v 
characteristic to the zero current axis. The half-gap 
value so obtained was corrected to absolute zero 
using the Muhlschlegel tabulation (the correction was 
negligible in this case). The result'Sgiven in the 
figure as AO),, and is also listed in Table 6.2, 
along with other half-gap values obtained in a similar 
manner. 

This extrapolation method for obtaining the 
energy gap value is the one prescribed by Giaever et al 


(1962). The prescription is somewhat arbitrary but 


experimentally convenient. 


B. Phonon Effects for Lead, Tin, Indium and Aluminum 


mg Lead. Phonon effects were first observed in 
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Pieure 5.1 
First derivative results for aluminum 


showing the gap edge singularity. 


43. 


sh 


bay | stuals “2 
_ ; 


iwitnls “ol eslugert avivevilteb je 7 
a 


«Yilislugole exgo6 gay edz 


% 
repre * 


—_ 


-20 


e) 


Al-42 (0:35°K) 
2 A0),= 420 LV 


20 40 60 
V-2A (Kv) 


80 


(HP@EO) Sp-ia 
Vay OSA «(GAS 


4, 


lead by Glaever et al (1962). These authors suggested 
that the variation which they observed in their first 
derivative characteristic was due to an energy-dependent 
Zap parameter, and they also pointed out that the 
variations occured in a voltage region near the Debye 
energy. 

Rowell et al (1963) subsequently published 
detailed first and second derivative experimental results 
for lead in conjunction with theoretical articles by 
Schrieffer et al (1963) and Scalapino and Anderson (1964). 

We have only studied two lead specimens, and in 
both cases the observed results were in excellent agreement 
with published results. The results obtained for one of 
these specimens are shown in Figures 5.2 and 5.3. 
Appreciable broadening of the gap edge singularity for 
lead is evident in Figure 5.2. The slight ‘overshoot' of 
the observed points relative to the BCS curve at 0.4 mv 
may be due to the superconducting probe distribution. 

The maximum amplitude value of 4.5 for the characteristic 
at the gap edge for this specimen is typical of the results 
obtained for most of the specimens studied here. The 
largest maximum value obtained for any specimen studied 

to date has been 20. These values are comparable to the 
few that have been reported in the literature. The 
broadening of the gap edge singularity may be due to the 


thin films used for specimens. 
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Figure 5.2 
First derivative results for lead 


showing the gap edge singularity. 
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The effect of the transverse and longitudinal 
phonon groups on the © characteristic for lead is 
readily apparent in Figure 5.3. The arrows correspond 
to predominant transverse and longitudinal phonon 
energies along symmetry directions for lead observed 
by Brockhouse et al (1962) with neutron scattering. 

The Lorentzian shaped curves shown below the experimental 
points are proportional to the phonon distributions used 
by Schrieffer et al (loc cit) in their model calculation. 
The structure at energies above 12 mv appears in 
Schrieffer's model calculation also and is therefore 
assumed to be sum and harmonic structure. 

e. Tin. First derivative results showing phonon 
effects for tin have not been published to our knowledge, but 
Rowell et al have published second derivative results for 
Sn-Sn0,-Sn tunnel junctions on two separate occasions 
(Rowell et al (1963) and Scalapino and Anderson (1964) ). 
Three tin specimens have been investigated here and the 
results obtained from two of them are shown in Figures 
5.4 to 5.7. Figure 5.4 shows the normalized first 
derivative results obtained from specimen Sn-40 at 2.2°K. 
Phonon deviations are clearly evident at voltages of 4 mv 
and 16 mv, and smaller deviations are discernible through- 
out the 4 to 16 mv region. The results shown in Figure 
5.5 were obtained with a different specimen and a 


superconducting probe distribution. It may be seen that 
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| ‘First derivative results for lead 


showing phonon structure. 
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at 242° Kvshowing phononitstrueture: 


Pieure Ss “5 
First derivative curve for tin 


at 0.35°K showing phonon structure. 
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the phonon effects are reproducible from one specimen 
to another, and that the superconducting probe 
distribution enhances the small variations. It is 
also clear that the phonon spectrum of tin is quite 
complex relative to that of lead, as one would expect 
from the crystal symmetries (white tin has an 
interlinked body centered tetragonal structure, while 
lead has a face centered cubic structure). The tin 
Spectrum appears to have a broad transverse group 

at about 4.5 mv, and a longitudinal group at about 

16 mv. These energies are 0.3W)p and 1.2Wpy 
respectively, where Wp is the Debye energy. 

Figure 5.6 shows the second derivative results 
obtained by taking first differences from the data shown 
in Figure 5.5. This method introduces undue scatter into 
the results, so that a visually estimated smooth curve has 
been Grawn in. The gross structure indicated by the 
arrows is in good agreement with the results published 
by Rowell et al (loc cit). 

A small portion of the original data obtained 
from specimen Sn-41 was re-analysed with a protractor- 
like instrument which enabled us to read slopes directly 
from the projected photographs. The slopes taken were 
estimated best fits over a 0.3 mv interval, so that 
noise present in the recordings with this period or 


less was in effect filtered out. The results were 


OH 


5 
if tet9so@ subnu 2eeouporluat bodyem aint .@s.2 surg th me 7 


Yitoeilb ascele beet ot eu Deldans ao LeMy 


ilosqe sno mon? sidigsbouges sis siosTie 
sdotq watvoubnoonsgva sit Jsdt bis _ 


‘. ,anolraivey Lf sme sds esonsanhe sok 


ar 
rm 
ci 


iosexs Bluow enc 86 ,bsel to gad? ot svideteg hail 
ym end mit stiaw) esintennyvse Cstagug ead aot? 
aftriw ,avudourte [snonsited betednss ybod bedatixeoad | 
ait sdT ,lerudtoutds oldys bersinss efa® @ aed Basle . 
quo Ts se rteyedaty bsowd 6 e¥sd oc Staedas muritooqe a, 
‘uocds ts quowyg Lantbutians! 6 ons y ve Gor Seage a 


gid 5.1 S48 pL) €.0 sus estguens SORRe wor 


pt eras 
«VEIenis, syed sdv eat gs s"tscdw _yievitoeqest - 


veaYT asyvizeriaeb bwoaSset soo ewoneg o.¢ 15g lt ‘  « 


te sisb stv moat asonensitib gatlt aniadaed ve bentstdo 
& 


avivo (Wwooms Sstkmigas yilfsualy «e Jans oe . etlueet 
ai? yo bedeolbnl stusstgvs seorg ed? «At wget need 
iefatiduq ativest odd dvi tnameetas boom AE ewes 


ty’ cat) £8 Pie acts! 


benisddo sJs6 Isntgiay sdt Io nofiaog Lisa 


-rodostiotg s diiw beeylaasear esw [Mune “a 


exeWw asisat seaols s”AT radqotgosoda re 
tad og ,fevissal Wm £.0 8 ‘Tevo set 4 
10 bolTs@.aidd dittw east: eases 
oxen eelueet si? sdhoD 


Figure 5.6 
Second derivative curve for 


tin - ao Qa 5a Kk. 


Pi gure ../( 
Expanded portion of Figure 5.6 


showingscritical. point curve fit. 
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plotted and formed a smooth curve. Tangential slopes were 
taken from this curve, and the results are shown in 
Figure 5.7 as points. A curve of the form —oa \Inlol 
has been fitted to these points (cf. page 1&8 ). The 
critical point energy obtained in this manner is 
Dia 2s. 050 mv "(1.3 Wp ) and it is in agreement with the 
value 17.7 = 0.1 mv reported by Scalapino and Anderson 
for a Sn=-SnOy-Sn tunnel junction. The absence of a 
well-developed singularity on the low energy side of the 
critical point (1i.e., xX negative) is also in accord 
with Scalapino and Anderson's results, indicating that 
the end point singularity is not sufficiently well 
isolated from other phonon effects. 

os Indium? some preliminary’ first derivative 
results showing phonon structure for indium have been 
published by Adler and Rogers (1963). Three indium 
Specimens have been investigated since then, and the 
results obtained from one of them are shown in Figures 
5.8 to 5.10. While the quality of the characteristics 
has varied slightly from one specimen to the next for a 
given type of specimen, the predominant phonon effects 
have always been found to be reproducible. The data 
for the critical point fit shown in Figure 5.10 were 
obtained in the same manner as the data shown in 


Figure 5.7. 
* Uncertainties quoted throughout this work are estimated 


outside limits for a given measurement. Insignificant 
figures will be denoted with a bar. 
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Figure 5.8 
First derivative results for indium 


at 0.35°K showing phonon structure. 
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Fagure 5.9 
Second derivative results for 


nam at 0, 35°R. 


Pigure 5.10 
Expanded portion of Figure 5.9 
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One may draw the following conclusions about the 
phonon spectrum of indium from Figures 5.8 to 5.10: 

a) the phonon spectrum is of a complexity 

intermediate to that of lead and tin, 

b) the phonon spectrum has a broad transverse 
group at about 5 mv ( 0.5 Wp) and a long- 
itudinal group at about 13 mv (1.4 Wp ), and 

c) the end point of the phonon spectrum has an 
energy of 14.8 £ 0.1 mv (1.56 Wp ). 

A critical point curve (Figure 5.10) has been 
fitted to both sides of the 'singularity' observed in the 
second derivative results at 14.8 mv. The fit to the low 
energy side of the singularity is somewhat uncertain, but 
the critical point energy coincides with that obtained 
from the curve fitted to the high energy side of the 
Singularity. It is thus possible to assign an approximate 
value of 0.3 to the ratio G/y (cf. pagei® ). The @ 
and & curves given by Scalapino and Anderson for lead 
Sugeest that this ratio is of the correct order of magnitude. 

4, Aluminum, Three specimens have been studied 
to voltages beyond the Debye energy, but phonon effects 
were only observed in one of these. The effects are 
very small and are only visible under ideal conditions. 
Figure 5.11 shows the first derivative results obtained 


for specimen Al-42*, The solid line labeled 'BCS' is 


* The author wishes to thank Mr. D. W. Penner for 
preparing and mounting this specimen. 
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Fagure 5.11 
First derivative results for 
aluminum at 0.35°K showing 


phonon effects. 
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the characteristic for two identical superconductors 

as calculated from the BCS theory by Adler (1963). 

The agreement between calculated and observed results 

is very good for applied voltages of less than 8 mv. 
Although it is perhaps speculative to base conclusions 
on the results obtained from only one specimen, the small 
deviations at higher energies are very probably phonon 
effects. The uncertainty in the © characteristic is 
somewhat less than the nominal * 0.0003 which we usually 
obtain with our present instrumentation. The excellent 
behaviour of the specimen is responsible for this low 
noise because it allowed us to use higher gain and the 
lower attendant bandwidth for recording. 

If one estimates the actual © characteristic 
from the observed data, then definite trends are in 
evidence at the arrows We, ,We,, and Wes ~. These 
energies correspond to the end points of the lower 
transverse, upper transverse, and longitudinal modes of 
the phonon spectrum for aluminum as obtained by Walker 
(1956) with diffuse x-ray scattering. At each of these 
energies the O& characteristic should be recovering 
from the minimum that occurs at a lower energy, and 
such a behaviour appears in Figure 5.11. Walker's 
results also indicate a strong peak in the phonon 
spectrum at an energy corresponding to the arrow 


labelled Wp . The @ characteristic is observed to 


= 2 


7 on 
Pe ir 
2thobAenvegque Isstinebt ows 107 otatredeneds ¢ 


: 


(£38f) YelbA ya ytoods 20d sdé movd boserualne a 
eifuzet bavisedo tris bstaivoiss neewied saomeeran | 


w 


ve 3 oss seo lo sopatiov betiiags tot boog vray er 
snotbuicoton sead of ovidstivoeqs adedteq at-st 
ifeme sdv ,nedlasda sto ylno mort benistdo eo tuner, edd 0 


noncig Yyldadorq ys) sx5 egigione tedgid 8° . 


et Stdafietoetsan DO edd at yonkerteon off 


vifeuep ew dotriw £000,¢ * fehioon wid nedd sacl 


$nsiiesxs 211] .UolgséAszastiuadteant taeszeta wo ies as 
iol 2iny so: sldienoges ei aenifoege of2 to wooded 
x4¢ Bos nt 54 tajfigid seu of su bewollie 22 saneoed eeton 
- 
metyioose "0% db iwingd Inannedss xwoh 
oitelivetestatio “D Lsutes ond astamites sao 7f ee 2 
it sts ebnet’ sticniteb aed? ,sisb bevissde — 
ses? . geld One , gold, go eworts odo ta 4 ca 
sewol pris lo esiteq bate sd ot Soogeeanes. 
© asbom lantbugtanc! one ,seTevensts Sega 
soniaW vo betitsido es muAtdiie ict muro eqe nonerg:. 
segs? lo dose SA .gnitedisge ysi-x oswllth Bgiw. sal | 


atitevpos: ad biyoite otjebrevosusdo S ade er sasagek 7 
has ,Vg tens tewol a 2 neon eae xe 


is = the _ 


e'teieW .f1.2 siygtl at evssqqs ayok 
femora af at pleas 


on i 


decrease rapidly as expected at this energy. 

Our results for aluminum are thus in good 
qualitative agreement with Walker's x-ray results. 
Tunnelling methods may provide more precise values 
of the end point energies than are available from 
x-ray methods, should a method of lowering the 
uncertainty in the observed OO results be devised. 
Modest improvements in the techniques of preparing 
Specimens and recording the results should accomplish 


the necessary improvement. 


C. Preliminary Results for Mercury 


It should be possible to observe phonon effects 


in tunnel junctions involving mercury, but such specimens 
are difficult to prepare since mercury melts at -40°C. 

One such specimen has been obtained in the manner described 
on page 35 but the results are unclear due to an 
excessive amount of filamentary conduction. 

The OO characteristic obtained for this specimen 
is shown in Figure 5.12. 

The violent flucuations in evidence are definitely 
not phonon effects, so that one may only conclude from 
Figure 5.12 that the phonon effects in mercury are less than 
0.03 in magnitude. This is not a very informative conclusion 
as we expect the phonon effects in mercury to have an 
amplitude of about 0.01 (see Table 5.1). However, the 


BCS trend is clearly in evidence, so that this method of 
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First derivative results for mercury 
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Specimen preparation should be capable of yielding good 


results. 


D. Semi-empirical Equations for the Strength and Temperature 
Dependence of Phonon Structure in Tunnelling Experiments 


Phonon structure experiments are not readily 
executed, so that it is desirable to know which metals 
might give phonon effects that can be resolved with a given 
piece of apparatus. The two limiting factors of primary 
importance are perhaps the lowest temperature that can be 
achieved in the cryostat and the resolving power of the 
apparatus used to obtain the dynamic conductance of the 
tunnel junction. Empirical knowledge at least of the 
strength of the phonon structure and its temperature 
dependence for a given metal in terms of available 
data would therefore be very useful. Reasonable 
available data might be taken as the Debye energy Op 
and the energy gap at absolute zero ey a . ht OYE? 
to set up an empirical equation for the strength of the 
phonon structure, consider the structure for lead. 
Schrieffer's model predicts a strong minimum in Oo." suet 
beyond CNT) ; Where Cp is the energy of the longitudinal 
phonon group. The exact nature of the minimum for 
Schrieffer's model depends on the ratio Wp/W-e , the 
electron-phonon interaction strength, and the Coulomb 
pseudo-potential. However, since Due also depends 


on the latter two quantities and since W2/ Ce will not 
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be widely different for various metals, one might expect the 
energy dependence of QO, and Li to be roughly universal 
from Schrieffer's model. Taking UW, = Wp and noting 
that 

Ti Pii+ at-Az 

ZU)* 
for 
Se Ae ae 


one might then expect 


Load 
where S, is the strength of the phonon structure at OK » 
Oia is the normalized dynamic conductance at the minimum 
of the predominant structure appearing at perhaps 1.2 Wy, , and 
C is a proportionality constant which we shall assume is the 
same for both strong and weak-coupling superconductors.* If 
one assumes further that the temperature dependence of A, and 
QO. is the same as that of the BCS half gap energy A(t) 


’ 


then 


Sea aeal (2ia)" 
ZS, Alo) (5.2) 


This equation gives the temperature dependence of S(t) in terms 


cee A eT) 


* The weak coupling superconductors are those for which 
No)V « 4 , where N&) and V- have the same meaning as in 
section II A. While the BCS theory is restricted to the weak 
coupling limit in some cases, Schrieffer's approach is 
believed to be free of this restriction. 
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of that given by the BCS theory for the energy gap, so that 
S&) may be exhibited directly as a function of the 


Se 
reduced temperature fd" eats 
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Since the energy gap in leadgutinsaindiumseand 
aluminum is developed essentially to its zero temperature 
Value at the lowest temperature available in the He3 
eryostat, it has been possible to evaluate C in equation 
5.1 for these four metals. The results are shown in 
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Observed Phonon Structure Strengths 


Metal S C 
( x 10-3) 
observed 
Pb 43 oo 
In 45 ta 
Sn Da 2 
Al Sas Sf 1 
estimated 
Hg ba5 2 


The value of So observed for an Al-Al specimen was 

Bieler. fone, This value is expected to be more than twice 
the value which one would obtain with an n-s system, and 
Since the values for lead, indium, and tin correspond to 
those for n-s systems, the value = UTRIO © has 
been entered for aluminum in Table 5.1. The values obtained 


for C may be seen to be essentially constant relative to the 


large variation of Se . 
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The value of Sg may be calculated for mercury 
using this value for C, and we see that the phonon 
structure in mercury should be readily resolved with 
the apparatus used here. 

The temperature dependence of S(t) has been 
observed for lead and indium. The results obtained 
are shown in Figure 5.14, and are in reasonable 
agreement with equation 5.e. 

Equations 5.1 and 5.2 would thus appear to be 
useful empirical guides for proposed phonon structure 


experiments on weak-coupling superconductors. 
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Figurei 5.13 
Observed temperature dependence 


of the first derivative curve for lead. 
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Experimental test of equation 5.2 
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CHAPTER VI 


TUNNELLING CHARACTERISTICS IN THE ENERGY GAP REGION 


A. Josephson Currents and Double Particle Tunnelling. 


Two mechanisms are known to exist for conduction 
within the energy gap (i.e., for V<A ). These are the 
Josephson current mechanisnl, which gives rise to a d.c. 
current for zero applied voltage and to a high frequency 
a.c. current for small applied voltages; and the double 
particle tunnelling mechanism which gives rise to a 
current jump followed by a conduction increase at each 
Of the voltages A, and Aj, . Both processes 
require that both metallic layers of the tunnel junction 
be superconducting, and are favored by a thin barrier 
layer. Since we have not used unusually thin barrier 
layers for the phonon experiments described here one 
would not expect either of these effects to be strongly 
evident in the gap region characteristics for our specimens. 
However, the current that appears in the gap region of 
the current-voltage characteristics shown in Figures 
6.1 and 6.2 may be attributable to double particle 
tunnelling if one assumes that the current jump at Da has 


not been resolved. The temperature-independent linear 


1. Josephson (1962) a 
2. Schrieffer and Wilkins (1962), Taylor and Burstein (1962), 


Adkins (1963). 
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Figures 6.1 to 6.4 


Current-voltage characteristics which 
exhibit double particle tunnelling 


effects. 
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increase in current beyond Ag, is in qualitative 


agreement with double particle tunnelling results obtained 
by Taylor and Burstein using a Sn-Tl tunnel junction. One 
may also see a further increase in conduction at 

V= Ameat+2 Aa, in Figures 6.1 to 6.4. This is the 
energetic threshold for a double particle tunnelling process 
in which two electron pairs are broken up in Mg and two 
excited pairs are produced in the aluminum, as indicated in 


Figure 6.5. 


Naar 2Da, 
Fermi levels 


Figure 6.5 


A possible double particle tunnelling process. The 
initial states are shown as open circles, and the 
final states as filled circles. Energy conservation 
requires that the energy gained by the two electrons 
in crossing the potential difference VY be equal to 


or greater than the energy of the final excited states. 
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It is not clear however why such a mechanism 
should produce effects in the current-voltage characteristic 
which are much stronger than the double particle tunnelling 


effects in the energy gap region. 


As one would expect, Josephson currents have not 


been observed. 
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B,. Filamentary Conduction 


Figures 6.6 and 6.7 show current-voltage 
characteristics in the gap region which display a current 
Jump at the origin as is expected for the d.c. Josephson 
current. However, the characteristics do not drop to 
zero just off the origin, so that some other conduction 
mechanism must be responsible. Giaever (1960) has 
explained a similar behaviour in an Al-A1,03-In tunnel 
junction by assuming a metallic bridge pierces the oxide 
layer. The bridge was assumed to contribute negligible 
conduction in the normal state. An order-of-magnitude 
estimate for one of our specimens based on the bulk 
resistivity of aluminum at room temperature gave a bridge 
Giameter value of 0.1 A. A bridge would only be conceivable 
if it had a diameter about two orders of magnitude greater 
than this, which is still plausible for this model. It 
does suggest that the bridges are very fine filaments. 

The behaviour of the Pb-39 characteristic shown 
in Figure 6.7 might then be explained as follows: 

The barrier layer was pierced by two aluminum 
filaments with different cross-sections. The resistance- 
current characteristic for the filaments is assumed to have 
5S broad transition at the critical current. As the voltage 
across the junction was increased the current in the 
filaments would remain relatively constant near the critical 


value, but once the filaments became normal the heat 
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Figure 6.6 
The current-voltage characteristic 
for a tin specimen showing filamentary 


conduction. 


Figure 6.7 
The current-voltage characteristic 
for a lead specimen showing filamentary 


conduction. 
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generated in them would prevent them from becoming 
Superconducting until the voltage was reduced well 

below that for which the filaments became normal. The 
Smaller of the two filaments in the characteristic shown 
in Figure 6.7 thus appears to have had an upper critical 
voltage of 0.2 mv, while the stronger one persisted in 

the superconducting state for voltages in excess of the 
gap voltage. It must eventually have become normal 
however since at large voltages the O&O characteristic 

for this specimen normalized to unity. The choice of 

two filaments for the explanation of Figure 6.7 is quite 
artificial. It would seem more reasonable to assume 

that there are either a great many filaments or none at 
all in a particular specimen. A multiply bridged specimen 
Should have a large number of small discontinuities through- 
Out its current-voltage characteristic. These small 
discontinuities should be apparent on our dv-V recordings 
because of the large gains employed for Sv sh Oo VERLS a 
in fact where we first encountered this type of specimen 
fault. Figures 6.8 to 6.10 show this behaviour for specimens 
Sn-36, Pb-39, and Hg-1l. The filaments in the tin and lead 
specimens appeared to be aluminum, since the current jumps 
at the origin of the current-voltage characteristics 
disappeared near the transition temperature for aluminum. 
Useable phonon characteristics were thus obtainable for 


these specimens at 2°K. The filaments in the mercury 


Ae 
oto 


=, 


gnimoosd mo) make dasverg bivow medt ak ; 

{Low besupe: saw egadlov edt Lida WD 

af? .Lemion sasoed atmemelit odd dobiw 40% 3 wet 
anwode otgatvedpsisdo edd al ednemeiti ows on% To° atten 7 

{sotvivo taqqu nB fer svet Of etsoqgs ewdd Jj. r) swath at hy 

nt betatesea eo segnette oct ofLtiiw vm S$. 0 to eustioy ~~ e 

eid to sesoxe ni seusifov toT sésse. beiiaemmnie sc 
[amron smoosd svsd yllavdneve teum JI .egetfLov 8 
giteinesdosiads © sts eogatfov syiat ts sonte tevewor 

to eotodo efT .vdtaw ot bextiamuion neomioegs etds 102 a a 

adiup at Y.2 etwalT Yo mottsasiqxs of? tol eomemall? ows ss) 


omuzes oF eidsnosset e1om mesa bivow 31 Matotitess : 


te sea 10 etnomeallt ysem isertn 6 t9ed¢le os eted? sedd 
nemtosge besbiad yiqttium A .nemtosge weilwotdiiaqg & at the 
-fyuorls astilunitnoos!b [lama to redmun sydel & eves biuvoda 
{isms sesiT .voitatresosisdo ogedLov-sneri9 att os . 
egnibtons: V-v8 two no dmetsqqs od bivoda setdtunt¢nooeth — a : 
aew etd? . VS tot beyolame antsy sag uel edt to _— . r 
nemtooiga to eqyt etds betesmoons Jaril ew st fost 
anamtoegs 101 tuotvaded eld’ wore 0L,9 of 8.3 wha ’ 
best bas ati eft al stmereilt oat -f-g3H bns ee 
aqme, tnmetivo edz sonta ooninats sd o¢ ferssiqh.iu 
gottetiedosieds egetiev-sastiwe odd to stab 3 38 
Atintowis tol e1uteareqmet moljtanst? edt vane 1qsetb 


7) aj ean auds eal sobiairosostsds 


Jan =i 


, 
ae st 
7 

7 


tes 


Figures 6.8 to 6.10 
dv-V recordings showing erratic 


behaviour due to filamentary conduction. 
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Figure 6:8 
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Figure 6-9 
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Figure 6°10 


Sn-36 
(0:8 °K) 


a] ad 


A 
BEE ok 
7] -_ 
| bb es 


4 * a ‘ 


8-2 ewplt | 7 


i- ph eS 
(A*SE-O) | 


edowe =, | 
5 : 


~-——e oo 
~~? ia aarp did 7 & 
i ! 7 : 


ec~-d9 
(A*3E-O) 


Or owe! 


: 8-H 
REO) 


136 


Specimen appeared to be mercury, and were quite strong, 


So that useable phonon characteristics were not obtained. 


oe Broadening of the Gap-edge Singularity 


The BCS theory predicts a singularity in the 
density of states at the gap edge of a Superconductor, 
and this should manifest itself as a jump discontinuity 
in the current at a. voltage V= BatAp fora 
tunnel junction consisting of two superconducting metals 

Ma and My - such a jump discontinuity is never 

observed, which indicates that the gap edge singularities 
in the actual density of states functions are 'smeared.' 
Giaever (1962) has investigated the gap edges using a 
Sn-Sn0,-Sn tunnel junetion, and has been able to get 
reasonable agreement between observed and calculated 
current-voltage characteristics when Hebel-Slichter 
smearing* is introduced into the BCS density of states. 

Other forms of smearing might produce agreement 
that is just ‘as good. In fact the major contribution of the 
Hebel-Schlichter smearing according to Giaever's calculations 
is to produce a linear current rise from zero instead of a 
jump at the gap edge and this result may be obtained simply 
by truncating the BCS distribution. The calculations with 
the truncated distribution are immensely simpler than 
Giaever's and we shall now outline them in order to compare 


the smearing observed in some of our Al-Al junctions to that 


ee SE ts Sts 


* See Adler (1963). 
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observed by Giaever for tin. 


In order to set up the truncation we denote the 


BCS distribution factor by G, where 


SOS ee 


Measuring energy from the gap edge in units of 


yields 
; Fr 
ep2 x) * , 
where X = Co= cs. 


The area under this curve between zero (the gap edge) and o 
; 
is SAW aC) 12 . Thus we are led to adopt the truncation 


form 


e) 5 Mm <-K& 
Ga (22) Pes a te 


(2x) 2%. oxeenaa . 


where Gy is to represent the smeared BCS density of states 
hactor.,. The current-voltage characteristic for a 
symmetrical tunnel junction* is then calculated as the 
convolution of Gg with itself. This produces a linear 
increase in the current when only the truncated portions 
OL Ga are effective, and the resulting normalized dynamic 


conductance is readily calculated to be 


C= os for 2A(iI-%)<V<2A . 


* A symmetrical tunnel junction is one in which the two metals 
are the same. 
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We can apply this result directly to our aluminum résults, 
Since we have always used aluminum as our base layer. Three 
Al-~Al specimens have so far been recorded in this laboratory. 
The current-voltage characteristic of one of them (Al-26) 

has already been reported by Adler (1963). Approximate 


values for od are listed in Table 6.1. 


Table ©, "7 


Values of the Gap Edge Broadening Parameter 


Specimen Reference 4 

Al=-Al 26 Adler 0.10 
Al-Al 34 This work 0.0 
Al-Al 42 4 0.03 
Al-400 A-Al-500 A Giaever (1961 0.02 
A1-400 A-A1-270 2 " 0.05 
Sn-Sn G55 
The @-values for Giaever's specimens were 


obtained by reading slopes directly from his published 
current-voltage characteristics. The dimensions shown are 
the fYinethicknesses;* "The eo -value which we have taken 
from Giaever's tin results agrees with the corresponding 
value for Hebel-Schlichter smearing which Giaever found to 
give the best fit to his experimental data. Thus our 
trucation parameter @& bears a reasonably close resembence 


to the Hebel-Schlichter smearing parameter 86 . The 
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rather large OC value for our Al-26 specimen may have 
been due to a.c. pick-up in our early measuring circuits. 
A.C. pick-up from the 60 cycle mains is always present 
and it causes our current-voltage characteristics to 
be averaged over a voltage width Vp about the applied 
voltage V , where Vp is the unwanted a.c. voltage. 
The peculiar grounding arrangement of our present 
circuitry has been evolved as a direct attempt to 
minimize Vp , which we believe is now about 1 pv 
or less. We should thus be able to detect ®& values down 
to 0.01 or perhaps less, but a problem will undoubtedly 
arise with the present instrumentation if gap smearing 
Le,t0 studied in detail. 

The mechanism for gap smearing is not clear. 
Giaever has suggested that strained films are the cause, 
while Adler suggests that the thinness of the films is 
responsible. Adler's suggestion would appear to receive 
some support from a recent theoretical article by Blatt 
and Thompson (1963). These authors find that the energy 
gap undergoes resonances as a function of film thickness. 
The resonance period is of the order of 10 A, and the 
gap averaged over a few periods becomes larger than the 
bulk value as the film thickness is decreased. Since 
our films are both thin and not uniform in thickness, we 
have here a possible explanation for both gap smearing and 


for the increased gap value observed for thin films by 
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Giaever et al (1962) and others. A direct check on the 
Blatt-Thompson theory might be effected by measurements 
on a series of tunnel junctions prepared with very thin 
metal films. A study along these lines is currently 


being undertaken by Mr. D. Penner in this laboratory. 


D. Some Further Energy Gap Values for Aluminum, Lead, 
gallo Indium, and Mercury 


Table 6.2 shows some energy gap values corrected 
to absolute zero obtained during the course of this work. 
Theeratl oO 2O0) , where T. is. the bulk transition 
temperature aya 2 SO) is the energy gap value at O°K, 
is also shown in this table. Some values reported by 
other workers for this ratio are included for comparison. 
It may be seen that the results obtained by various workers 
are not consistent to within the stated uncertainties in 
measurement. Part of this inconsistency may be due to the 
use of thin film specimens, since Giaever et al (1962) 
and others have found that a decrease in film thickness 
causes an increase in the energy gap and the film transition 
temperature. A more suitable transition temperature for 
the ratio £90) might therefore be the value obtained 
directly from the tunnelling characteristics*. However, 


Douglas and Meservey (1962) have found the ratio 240) 


Cc 


ee 

* The transition temperature determined from the resistance 
of the film itself is probably meaningless due to edge 
effects. 
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Table 6.2 


Energy Gap Values 


Reference Metal 240) 22° 
(mv ) i a ae 

Giaever et al. (1962) Pb 2,08 4.33 

Sn Tn) mic) 

In 1.05 ayty 

Al 0.40 SeoT 
Douglass and Meservey (1962) Pb i 4.20 

Al - 3,3 * 
Zavaritskii (1961) Pb 2.66 4.26 

Sn ieee. 3.47 

In igce 3.45 

Al - ie A 
Bermon and Ginsberg (1963) He 1 abdast.0.04 4.5 
This work Pb Ac Bie. Oe Nee) 

He 1.63°5-0708 4.51 

Sn 170 0.05 <61. 


3 
In eOSa eds 3.75 
\ 


Al i ae = 002 


* Calculated on the basis of the film transition temperature 
which varied from 1.3 to 1.5°K for different specimens. 
The remaining ratios are based on the bulk transition 
temperatures. 
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where Ye is the "film transition temperature", is still 
dependent upon film thickness for aluminum films. We have 
accordingly used the bulk transition temperature for our 
faG.os, shown in Table 6,2. 

This convention which is also adopted by Giaever 
et, al., is experimentally convenient, since temperatures 
need not then be measured with precision. 

The ratio 2.06) given in Table 6.2 for 
aluminum is abnormally high. In this case the transition 
temperature of the films as observed from the tunnelling 
eneracteristics was 1.40 © 0.05°K, giving a ratio 

2660) = 3.5404 . This rather large discrepancy 
a ae the film transition temperature and the bulk 
transition temperature for aluminum appears from the 
literature to be a characteristic of that metal in 


that the corresponding effect is observed to be far 


less pronounced for lead, tin, and indium. 
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CHAPTER VII 


CONCLUSION 


A, Results for Lead, Tin, and Indium 

The phonon structure effects observed in lead and 
tin have been in agreement with similar results published 
in the ldterature, while our results for irtdium cannot 
be compared with other published data, observed or 
calculated, in that we have not been able to locate any 
such data. Some end point energies for the phonon spectrum 
of aluminum have been estimated from our aluminum results, 
and our estimates agree with Walker's x-ray results to 
within the uncertainty of about 3% in the x-ray data. 

The values which we have obtained for broadening 
of the gap edge singularity and for the energy gaps of the 
various metals studied are in good agreement with the 
corresponding values which have appeared in the literature, 
indicating that the specimens which we have prepared are 


comparable to those prepared by other workers in this field. 


B. Thermometry Results 


Some preliminary work has been done in the use of 
tunnel junctions as thermometers, and the results appear 
quite promising in that the observed discrepancies are 
believed to be due to the experimental arrangement, and 
not the tunnel junction itself. For example we find that 


it is not possible to lower the temperature of the He? 
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chamber below 0.5°K when the indicated vacuum space 
Picoeure 18 above 5 x 107-6 mm Hg., but that temperatures 
of 0.35°K may be attained for a pressure of 1 x 1076 mm 
Hg. Since the specimen chamber has a surface area 
comparable to that of the He chamber, it is quite 
plausible that there is still enough heat leak to the 
specimen chamber at an indicated pressure of 1 x 10-6 


mm Hg. to cause our observed discrepancies in 


temperature. 
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APPENDIX I 
APPARATUS FOR MEASURING CHARACTERISTICS OF SUPERCONDUCTING 


TUNNEL JUNCTIONS 


The notation 4wv_ has been used in this paper 
for one of the voltages observed. This voltage has been 
referred to as Sw _ in the preceding text in order to 


avoid possible confusion with the gap parameter A 
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Apparatus for Measuring Characteristics of Superconducting Tunnel Junctions* 


J. S. Rocrers, J. G. ApLER, AND S. B. Woops 
Department of Physics, University of Alberta, Edmonton, Alberta, Canada 
(Received 8 October 1963; and in final form, 4 November 1963) 


Much information about the density of electron states in a superconductor may be obtained using electron tunnel- 
ing techniques. It is of particular interest to measure the normalized dynamic conductance of superconducting 
tunnel junctions as a function of applied voltage. This quantity is the ratio of the dynamic conductance (di/dv), 
when one or both metallic members of the tunnel junction are in the superconducting state, to the dynamic con- 
ductance (di/dv), when both members are normal. A new method has been developed which enables measurements 
to be made of (di/dv)./(di/dv), to a few parts in ten thousand. With this method only a bridge circuit, a galvanom- 
eter amplifier, and an oscilloscope are used. The galvanometer amplifier has a passband from de to a few cps and 
an input noise voltage of about 3X10~* rms V. The circuits of the bridge and amplifier are presented and analyzed. 
The operation of the circuits for measuring the characteristics of low resistance specimens (S500Q) is described. 


INTRODUCTION 


UCH information about the density of electron 

states in a superconductor may be obtained using 
lectron tunneling techniques.!~* It is of particular interest 
0 measure the normalized dynamic conductance of super- 
conducting tunnel junctions as a function of applied 
oltage. This quantity is the ratio of the dynamic con- 
uctance (di/dv), when one or both metallic members of 
he tunnel junction are in the superconducting state, to 
he dynamic conductance (di/dv), when both members 
re normal. Those variations in the effective density of 
lectron states in a superconductor which mirror the 
honon spectrum can be extremely small in some metals,4 
ut they have been described theoretically? and can be 
harted experimentally with sensitive conductance meas- 
rements. Two methods of making conductance measure- 
ents with the requisite sensitivity have been described!” 
cently. The bridge method described here, although 
mewhat less direct, is much simpler to construct and is 
pable of yielding the normalized dynamic conductances 


*This work was supported by the National Research Council of 
anada. 

t Present address: Department of Physics, Dalhousie University, 
alifax, Nova Scotia, Canada. 

mo H. R. Hart, and K. Megerle, Phys. Rev. 126, 941 
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tters 10, 334 (1963). 
*J. R. Schrieffer, D. J. Scalapino, and J. W. Wilkins, Phys. Rev. 
tters 10, 336 (1963). 
‘J. G. Adler and J. S. Rogers, Phys. Rev. Letters 10, 217 (1963). 
*D. E. Thomas and J. M. Klein, Rev. Sci. Instr. 34, 920 (1963). 
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of suitable specimens with an accuracy of a few parts in 
ten thousand. 

A typical tunnel junction between two normal metals 
is not strongly nonlinear in the region of interest, since the 
junction resistance only decreases by a factor of two as the 
applied voltage is increased from 0 to about 250 mV. If such 
junctions are cooled so that one of the metals becomes super- 
conducting, a strong nonlinearity develops at the origin 
in the 7-v characteristic due to the energy gap in the density 
of electron states of the superconductor. Also, far weaker 
deviations from the characteristic obtained for normal 
metals occur at voltages corresponding to the energies of 
the transverse and longitudinal phonons in the metal 
lattice. It is these small deviations that are of specific in- 
terest here. They generally occur in the voltage range be- 
tween 5 and 50 mV where the 7-v characteristic is approxi- 
mately linear. This near linearity is exploited by our 
bridge but, since highly nonlinear regions can be approxi- 
mated as closely as one wishes by a succession of sufficiently 
short linear relations, the apparatus may be used, with a 
little extra patience, to investigate the energy gap region 


as well. 
THE BRIDGE 
Principle of the Method 


In Fig. 1 a bridge is shown, one arm of which contains 
a nonlinear passive element R,, having a voltage v across it, 
As long as the input resistance R, of the preamplifier is 
large enough, the current through it will be negligible and 
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determination of the current i through R, becomes a 
load-line problem involving the i-v characteristic of R,, 
the bridge voltage V, and R,. The right-hand-side of the 
bridge may be viewed as a load-line problem as well; in 
particular, if the two upper arms of the bridge are equal, 
the two load lines will coincide as shown in Fig. 2. If, in 
addition, a constant voltage V, is subtracted from the 
bridge unbalance voltage, a new voltage Av is obtained. A 
suitable choice of Ra and V, will permit Av to be the voltage 
difference between the 7-v characteristic of R, and a tang- 
ent line to it, as indicated in Fig. 2. It is clear that Av will 
be small in the region of vp and so may be amplified con- 
siderably in presentation. A sequence of such Av—v pres- 
entations thus becomes a sensitive map of the nonlinearity 
of the specimen R;. 


Bridge Equation for a Three Terminal Specimen 


The bridge in Fig. 1 uses a three-terminal connection 
to R;. In practice, however, four terminal connections to 
R, are necessary as lead resistance changes can be appreci- 
able when the specimen changes from the normal to the 
superconducting state. However, for simplicity the analy- 
sis of the three-terminal system will be considered first. 


Ra weer 
Fic. 1. Bridge circuit jee NT 
with 3-terminal connec- a Woy Se 
tion to Rx, the non- Ec poh 1389 
linear resistance. TT cower ice ws 
Ea" - Preamplifier 


Let R, in Fig. 1 include the resistance of the upper cur- 
rent lead to R,, and assume R, to be large enough that the 
resistance of the upper potential lead to R, may be 
neglected. On the lower end of R, only one lead is used 
whose resistance then is part of R,. On writing three suit- 
able loop equations and solving for 7, one obtains 


t= R(RoRa)'{v+ [1+ RaR1+ (Rat R)Ra(R,R) Av} 
+(R+Ra)(RaRa) Vy. (1) 
This expression may readily be differentiated with respect 
to v. If we append subscripts n and s to distinguish be- 
tween measurements on R, when both metals are in the 
normal state and when one is in the superconducting state, 
respectively, then for Ra<R, and Ra~R, we may write 
for the normalized dynamic conductance, 
(di/dv), 
(di/dv)n 
Uscpicaat Hee [1+ (Ras/R)+2 (Ras/Ry) |[d (Avs) /do ]} 
RosRas{1+[1+ (Ren/R)+2(Ran/ Ry) Id (Ary)/dv]} 
(2) 


Fic. 2. Superimposed 
load line graphs. 


For regions not too close to the energy gap Ran Ras an 
if R,>10Ra, we can expand (2) and obtain 


g~1+[d(Av,—Av,)/dv | 
+[ (Ran—Ras)/ Ras |-+ [Ran Ras)/ Kaa © 


Let R,, be the resistance of the upper current lead when tk 
specimen is in the normal state. Then it can be seen frot 
(3) that the error 6,,g, introduced in g due to R,, becomin 
zero (superconducting) will be 6,g=R,/R. This error ca 
be made negligible by choosing a large enough value for 
Similarly, denoting the lower current lead resistance b 
Rj, then if it becomes superconducting, Ran— Ras will 
in error by an amount R; so that 6;g~R1/Ran. This errc 
can not be made negligible in practice since there is a 
upper bound to Ran imposed by leakage conduction acro 
the specimen and by requirements of the preamplifier t 
be described later. 


Bridge Modification for a Four Terminal Specimer 


As indicated above, difficulty is encountered in tk 
bridge of Fig. 1 when part of R; becomes superconductin; 
This problem is avoided when R; is the normal lead of 
specimen with only one metal superconducting, but mea 
urements in which both metals are superconducting al 
also of interest. The bridge shown in Fig. 3. avoids th 
difficulty by making R; effectively part of Ry. In order 
see how resistances of the potential leads affect the bridg 
operation, let us first consider a bridge which has bee 
balanced with R;=R,=0. This balance condition, 
course, is dependent on the current through the nonline< 
specimen resistance. If R; and R, become nonzero, a smé 
unbalance current will flow through the lower potenti 
lead R,. Neglecting changes in the other loop currents, or 


Fic. 3. Bridge cir- 
cult with four-ter- 
minal connection to 
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‘is eter preamplifier. 


Galvanometer 


' 


obtains from the lower current loop 


his unbalance current causes an error in the measure- 


ment of Av of R,67, which is equivalent to a change in 
Ra of 
Oh R,61/i~ RiR_/2R, 


ko that the error in g due to a change in R; is now 


; 5.g=R,6Ri/2RRz. 


if 

| his error may be made negligible. 
# It is essential for the Av-loop indicated by heavy lines in 
Fig. 3 to be free of thermal voltages® since the value of Av 
“Ms required with a precision of about 10~-* V. Thermal 
Woltages elsewhere in the bridge do not affect the measure- 
@mnents because the specimen resistance and Rg are small 
‘Fompared to the other bridge resistances. It should also 
be remarked that the low impedance of the Av-loop is 
#ssential to the operation of the preamplifier described in 
he following section. For this reason a Kelvin bridge is not 
suitable for these measurements. 


THE PREAMPLIFIER 


From the foregoing it is clear that the amplifier must 
ave a low noise level, and the input impedance should be 
igh. The passband must include de but need not exceed 
n few cycles per second. A modification of the galvanometer 
hmplifier described by MacDonald’ has been developed 
hich satisfies these requirements quite well. A simplified 
Hiagram of this amplifier is shown in Fig. 4. 

The optical arrangement shown differs from the sym- 
etric arrangement recommended by Jones® only in the 
placement of the photocells. Here the photocells are placed 
front of the symmetry position of the lamp. Lens Li 
orms an image of the lamp filament on the galvanometer 
mirror. Lens L2 forms an image of mask M; in the plane of 
nask M». This image is a vertical rectangle of light a little 
ider than Mp» and positioned so that some light passes 
ach side of M». Since the lens arrangement is symmetric, 
3 would cause the light passing M2 to reform into an 


SF. K. Harris, Electrical Measurements (John Wiley & Sons, Inc., 
ew York, 1952), pp. 180-185. : 
™D. K. C. MacDonald, in Handbuch der Physik, edited by S. 
liigge (Springer-Verlag, Berlin, 1956), Vol. 14, p. 159. - 

#R. V. Jones and J. C. S. Richards, J. Sci. Instr. 36, 90 (1959). 


» ADLER, AND WOODS 210 


image of the lamp filament at the symmetry position of 
the lamp were the light beam not intercepted by the 
photocells. The position of the photocells is chosen so that 
two small, separated, vertical strips of light fall on them. 
The widths of these strips are modulated by the motion of 
the galvanometer mirror. Since the photocells are con- 
nected in parallel but with opposite polarity, they generate 
either a positive or a negative signal which is fed back to 
the input after amplification. 


Analysis of the Preamplifier Circuit 


(a) Input impedance. Since the galvanometer acts as an 
error sensing eleme.t, and since 1;<t2, one obtains for the 
de input resistance, 


W ie = Av/14 tod Rpi2/y= aR ,, 


where a is the current gain. Typically, a~104 and R; is 
between 20 and 2000 Q, so that R, is of the order of 10° . 

(b) Gain. Retaining the approximation that negligible 
error signal appears across the galvanometer, one obtains 
for the over-all de gain 


A = €o/Av= (Rot Ry)/R;. 


Since Ro= 200 kQ, the gain is very nearly linear in R;. 

(c) Dynamic stability, vibrational noise. A good deal of 
information may be obtained about the galvanometer 
amplifier by determining its transfer functions in the 
open-loop condition.’ One may show from the differential 
equations for the closed-loop circuit of Fig. 4 that the 
currents 7; and 72 are summed in effect in the resistor Ry, 
so that the equivalent open-loop circuit is that shown in 
Fig. 5. The corresponding block diagram for the closed- 
loop system is shown in Fig. 6. In this diagram e, is a 
voltage equivalent of the noise induced by vibrations and 
introduced into the amplifier by horizontal rotations of 
the chassis relative to an inertial system. It may be evalu- 
ated from the differential equations for the closed-loop 
system. If / is the moment of inertia of the galvanometer 
coil, Rin is the total resistance of all elements in the 7;-loop 
indicated in Fig. 5, ¢ is the flux linkage of the galvanometer 
coil, D is the differential operator with respect to time, 
and @» is the angle the amplifier chassis makes with a suit- 


5 j Fic. 5. Open-loop circuit for the preamplifier 
R¢ Ss R¢ $ Re Cc 
rae es { 


ff oe 


9J, J. D’Azzo and C. H. Houpis, Control System Analysis and 
Synthesis (McGraw-Hill Book Company, Inc., New York, 1960). 


211 SUPER CON DUGD ENG YL NONE Lae JUAN Crom 


6, Bo 
—+| Go 


AV se ‘a 
Fic. 6. Block dia- rep ne Olt 


gram for the closed- at 


loop system. Hs | R¢ ig 
3 n 


able inertial frame of reference, then 
€n= RinI DO2 ‘p= RinKDO2/w17¢, (4) 


where w, is the natural resonance frequency of the gal- 
vanometer in rad/sec, and K is the stiffness constant of 
its coil suspension. It is clear that if Rinx< Rea, where Rea is 
the damping resistance specified for the galvanometer, 
then G;, is the transfer function of a highly overdamped 
second-order system. Denoting the damping constant for 
the G; system by &, one obtains 


Git = (KRin/ e) A+ 2EDon t+ D*an™) 
= (KRin/p) (1+ 2&Dwr)[1+ (26) Don], (5) 


where 2 V2R.a/Rin>1. Since R-a=10 kQ for the gal- 
vanometer used here, the overdamped condition is appro- 
priate. If the frequency dependent loading of the RoC 
network on the photocell post-amplifier A; is neglected, 
a suitable choice of 6;, the angle that measures the rotation 
of the galvanometer mirror relative to the amplifier 
chassis, may be made so that e9= 61, where & is a constant. 
Thus Go=k. 

The forward transfer function of the system is then con- 
stant to a corner frequency w; with the value 


Wi=Wn (2) Oi1@, Rains Rea, 


beyond which it decreases with a slope of 20 dB/decade to 
ws, the second corner frequency, located at 


Wo= DEO ge 1.40,Rca/Rin- 


Beyond w». the transfer function decreases at a rate of 
40 ‘dB/decade, with an attendant phase lag of 180°. Con- 
sequently, phase lead compensation must be introduced if 
the closed loop system is to be stable. The network shown 
in the is-loop of Fig. 5 is a standard phase lead compensa- 
tor, so that 


H=B8(1+7D)/(1+6rD), 


where B=R;/(Ro+R,) and r= R2C. Lead compensation is 
introduced for frequencies nominally between w,=7~! and 
wy= (Br), but is effective over a much wider range than 
this for the small values of 6 obtained here. 

The gain of the closed-loop system is 


Ax~1/H, (6) 


and has the value of 8! up to the corner frequency wa 
after which it decreases at a rate of 20 dB/decade. The 
actual gain characteristic is somewhat more complicated 
than Eq. (6) indicates because of the various approxima- 
tions made in the derivation, but it is a satisfactory guide 


to the amplifier performance. The approximations au 
excellent for the de behavior. 

(d) Noise level. Vibration appears to make the Haid 
contribution to the noise level of the amplifier. From th 
central term of Eq. (4) it will be noticed that the vibratio) 
noise voltage e, is independent of the stiffness constant ¢ 
the galvanometer suspension, so that using a rigid sus 
pension does not help. From the right hand side of Eq. (4 
it is evident that a galvanometer should be chosen fo 
which M=7°R;,/S is a minimum. In this merit figure, 4 
is the free period of the galvanometer and S is its curren 
sensitivity expressed in, say, mm/yA at one meter. Th) 
use of a high resistance galvanometer increases the toler! 
ance of the merit figure to high resistance specimens. 

It is also evident from Eq. (5) that a galvanometer witl 
a high current sensitivity is desirable so that the photocell! 


negligible when referred back to the input. 


OTHER FEATURES OF THE CURVE TRACER 
Circuitry 


The detailed circuit shown in Fig. 7 contains some 
further features that require comment. Although standar¢c 
shielding and construction practice were followed, som 
undesirable parasitic effects remained until the balance 
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Fic. 7. Bridge and preamplifier circuits. Heavy lines denote thermal- 
free wiring. Unless otherwise specified, all resistors are  W, 5% toler- 
ance, all resistance values are in ohms, and all capacitance valucs 
are in uF. 
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Fic. 8. Bridge supply circuit. Ri, Ro may be adjusted stepwise to 
e following values: R; =0, 470; R.=0, 390, 1.5 K, 4.7 K, 15 K, 47 K, 
0 K, 470 K. The electrolytic capacitor value is in uF /dewv. 


ntrol and the symmetrical configuration of the pre- 
plifier were introduced. 

A Tektronix 502 oscilloscope has been used for recording 
e results. This oscilloscope permits differential dc inputs 
both X and Y amplifiers, and has a resistance of 1 MQ 
nnected from each input terminal to ground. The effect 
these resistances on the bridge operation, particularly 
the grounding of undesired signals, must be considered 
en choosing the input modes to use. We have used 
fferential inputs to both channels of the oscilloscope and 
en provided an ac ground for the bridge through the 
pacitor associated with the balance control. Without 
is a 60 cps signal that was not entirely rejected by the 
cilloscope input appeared on the entire bridge circuit. 
e ac ground is adjusted with the balance control unit 
til it is at the virtual dc ground point opposite the 
Ecilloscope connection to the bridge. The Av-v trace is 
ependent of the sweep rate for slow speeds” when this 
justment is properly made. It is readily seen that with 
is arrangement feedback may occur from the output of 
e preamplifier to its input via the oscilloscope input re- 
tance, the balance control, and the bridge circuit. This 
dback is minimized by the symmetry of the circuit 
sign. 

It is interesting to notice that a complete symmetriza- 
n of the preamplifier would somehow involve a second 
lvanometer. Another advantage might be gained with a 
ond galvanometer by arranging it so that vibration 
ise would contribute to a common mode signal and thus 
rejected. A mechanical filter!’ has been tried to reduce 
effect of vibration, but it is difficult to eliminate long 
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IK 
v Rs Vv Fic. 9. Vy potentiometer circuit. 
el ion out R; may be adjusted stepwise to 
ru Ov22; 15; or 150 K. 
100 10 turn 
Vy vernier 


Slow speeds are those for which the preamplifier can follow the 
racteristic without lag. 

' The preamplifier chassis has been suspended from a spring, so that 
4s a free period of about 1 sec for vertical oscillations, and of 
ut 5 sec for rotations about the vertical axis. Oil damping is 
vided for both these degrees of freedom, and a shroud covers the 
Ssis and suspension to keep out air currents. 


. 


amy te 


ROGERS, ADLER, ANDXWOODS 


bo 
— 
bt 


v (mv) 


I'te . 10, Av—v traces obtained from the oscilloscope. 


period noise in this manner. The 30 cps “twin-T” notch 
filter in the output of the preamplifier suppresses a vi- 
brational resonance associated with the galvanometer 
suspension. 

The purpose of the transistor amplifier is to extend the 
frequency range of the preamplifier. The transistors are 
operated near cutoff (zero bias) so they require very little 
power. Some drift results from this unconventional biasing, 
but it is negligible compared to the photocell drift. The 
zero adjustment in the transistor amplifier has only a 
slight effect on the preamplifier zero because of the large 
gain which precedes the transistors, but it does have a large 
effect on any common mode ac signal from the photocells. 


Operation 


Sweep. The oscillator is used mainly for monitoring the 
specimen characteristics between measurements. The use 
of the manual sweep instead of the oscillator frees the cir- 
cuitry from a mains connection, and it is then electrically 
much quieter. The position (range and vernier, Fig. 8) con- 
trols are used to set the high voltage end of the sweep 
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Fic. 11. Normalized dynamic conductance results obtained from 
the traces shown in Fig. 10 (open circles) along with results obtained 
from similar adjacent traces (filled circles). The experimental scatter 
is evident in the overlap region. The characteristic approaches unity 
to within this experimental error at higher voltages. 
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range and the voltage is swept downward from this value 
with the sweep vernier. The sweep range is adjustable in 
half-decade steps so that any desired sweep may be made 
to occupy three or more turns of the vernier potentiometer. 

Av-v tracing. The preamplifier must not be subjected to 
large voltages since the galvanometer receives the full in- 
put voltage if the light beam is driven off the photocells. 
To maintain the preamplifier output within reasonable 
bounds the Raz and V, controls are adjusted in conjunction 
(Fig. 9). Variation of Ra causes the Av-v pattern on the 
oscilloscope to skew and shift vertically; V, is adjusted 
to oppose the vertical shift. In Fig. 10 a typical Av-v plot 
is shown and the analyzed results obtained from it appear 
in Fig. 11. The plot is simply a photograph of the oscillo- 
scope face. The information is retrieved from the photo- 
graph by comparison with a raster generated on the same 
oscilloscope, so that distortions in the various optical 


systems utilized cancel out. The data from the Av-v p] 
are processed by a digital computer according to Eq. (/ 
Little resolution is lost by recording on the oscillosco 
screen since for a typical preamplifier gain of 10%, a uni 
change in the slope of the pattern only causes a change 
the normalized dynamic conductance of 10~* assuming th 
equal gains are used on the vertical and horizontal osc 
Joscope channels. 

i-v tracing. For cases where negative resistance regio! 
are not required, the galvanometer amplifier may be us¢ 
to preamplify v while a large value of R, yields a lar; 
v;. If a negative resistance region is to be traced, the voltag 
supply for the specimen must be of low impedance. Or 
thus chooses a value for R, smaller than R, and uses tl 
low-z mode of the manual sweep. The galvanometi 
amplifier may then be placed across R, to amplify tk 
small v; which results. 
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APPENDIX II 
BRIDGE CIRCUIT FOR DETERMINING THE CAPACITANCE OF TUNNEL 


JUNCTIONS 


An extended Wheatstone bridge with a phase sensitive 
error signal display has been found useful during the 
preparation of tunnel junctions. 


Lteebasic circuit is shown in Figure A 2.1. 


Figure A 2.l 


Basic bridge circuit for capacitance measurements. 


fects tanure, C=C. exp. gut is a 
sinusoidal voltage applied to the bridge CL Vrculy,, Gi es. 
is the error signal, and the combination R,- Cy is to 
represent the specimen. One may readily derive the equation 


GQ. a, 1+u*t (tasaar) + EOtieyen , (4 2.2) 
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where on ae 
\ ae 
(eg = eee } 


C2 = Rez ) 
ar = ‘Gy =i 3 
and G is the frequency in radians per second. 
It therefore follows that the null conditions for 


the bridge are 
B= a and AT=O , 


or 
equivalently, LA 252) 
Re a ee 
32 
and oy te 3 2. ort (Ra 24) 
\-B2 


where (3 = a § 


While these null conditions may be more readily 
obtained by conventional methods, the derivation given above 
is more closely allied with the actual measurement operation 
which is now to be described. 

Figure A 2.2 shows the operational form of the 
bridge circuit. It may be seen that the horizontal 

deflection of the oscilloscope display is proportional 
to e » while the vertical deflection is proportional to 
Seg . The displayed pattern is therefore a straight 
line for low frequencies ( WtTKA ) as may be seen 


from equation A 2.2. Suitable values of Ro and (2 will 
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Figure A 2.2 


A bridge circuit for measuring the capacity of 
a tunnel junction. Ré¢ may be adjusted step- 
wise to the following values in ohms: 10, 31.6, 
LDOKGLHS LK SS LOK LOK eter ORs IOUK SU UR 
1M,0o0 . The resistors used for Rg are 
1%, 1/2W carbon film resistors. The 
capacitor Cs is a Heathkit DC-l, and may be 


adjusted from zero to 0.1mMf in 100 Bef steps. 


7 
% 


R3z: 100 ohm, 10 turn, 1/4% lin. 


% 3 Hewlet- Packard 200 CD 
Oscillator 
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make the straight line horizontal, and the null condition 
A 2.3 is then satisfied. An increase in frequency will 
cause the displayed pattern to open into an ellipse as a 
result of the term gw Er oe io Oi in equation A2.2, 
so that the value of C2 which causes the ellipse to 
close again satisfies the null condition A 2.4. This 
method of null attainment is suitable for specimens having 
a characteristic frequency W = ee which is at the high 
end of the operational frequency range of the bridge 
(66 hp. See+0)150,K,C.5. for the circuit shown in Figure A 2.2.) 

If the specimen has a characteristic frequency which 
is at the low end of the operational frequency range of the 
bridge, ait. will,»not,be,possible to,.satisfye ther condition 
Gi tec<el wedn this case,thesinitwal display» pattern 
will be a tilted ellipse for the lowest frequency available. 
If one then increases the frequency to 50 K¢-S- ,the ise 
terms in equation A 2.2 should be dominant... Suitable 
adjustment of Rea and © will cause the displayed pattern 
teecose, intosa.strateht.dine,rand finally Co may be adjusted 
to make the straight line horizontal. The null equations 
A.2.3.and.A 2.4.may,.then:berapplied as. before. »It.will, be 
noticed that the null conditions are independent of frequency, 
so that once a null is obtained, no matter by what means, it 
should be independent of the frequency. 

The above analysis is based on the ideal circuit 


shown in Figure A 2.1, but the actual circuit is a far more 
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complicated situation which behaves only approximately as 
the ideal circuit. We now consider the steps taken to 
bring the actual circuit behaviour into line with the ideal 
behaviour. 

It will first be noticed that an attenuated 
Slgnal is fed to the horizontal amplifier of the 
oscilloscope. The attenuation is chosen so that both 
channels of the oscilloscope may be operated at the same 
gain, then there is little relative phase shift introduced 
by the oscilloscope. 

A shielded lead, having an appreciable capacitance 
that must be added to the readings of the decade capacitor 
box, connects the capacitor Cg to the rest of the circuit. 
Finally we notice that the resistor Ra is of the helical 
Coll type, sa that it has an appreciable distributed 
inductance. It has been found convenient to correct for 
this inductance by making an additional correction to 
the value Cy . The net correction ac CG.d ie: has 
been obtained as a function of C, , Rg , and @, , where 

Rs ds themvalue of Re corrected for the input 
resistance to the ground of the 'A' terminal of the vertical 
channel of the oscilloscope. (A suitable dial on the 
resistor allows @2 to be read directly to £ 0.14). 

The correction dc has been obtained simply by inserting 
known values of C, and R, into the circuit. 


A typical set of correction curves for ee is 
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given in Figure A 2.3. Two apparently identical bridge 
constructions have yielded somewhat different correction 


curves, so these curves must be obtained for each bridge 


circuit made. It is implicit in the nature of the correction 


that OC should be small compared to C, . Our 
typical specimens have C, = 30 m pt » and dC 
is usually about 2mpf  teeyvy bout Ss of.) 


If the values of 5c are appreciable compared to C, 
the specimen capacitance cannot be determined with this 


DrLigescircuit, 


Rs rik, OC=+ 0-22 
316 
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Figure A2e.3 


Capacitance corrections. 
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APPENDIX III 


TUNNEL JUNCTIONS AS THERMOMETERS 


The current-voltage characteristics of a tunnel 
Junction in which at least one of the metallic films is 
Superconducting always possess information relating 
directly to the temperature of the Superconducting member 
(or members), so that such tunnel junctions may be used as 
thermometers. 
As a simplified example of this application, let 
Ra) = Ite be the density of states distributions for a 
Symmetric tunnel junction. The distribution diagram 
corresponding to the tunnel current integral for an applied 
voltage slightly less than the energy gap voltage* is shown 
in Figure A 3.1, along with the i-v characteristic in the 
gap region for this same tunnel junction. The current 
i(2a°) is due to thermally excited electrons in the 
distribution AN tunnelling into the distribution Pe, 
and also to the thermally excited holes in the distribution 
Ps, tunnelling into the distribution Bo. « Using 
the approximation of a constant tunnelling probablility, one 


may write an approximate expression for (toh) directly: 


; Pate y 
Peery ec fe * (&)* 4-06 des 
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The first term in the integrand is the Fermi function 
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Figure A 3.1 


Illustration of the use of a tunnel 


junction as a thermometer. 
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approximation for + >» 1 


the asymptotic form of the BCS function G for W=aA* ; 


» the second term is 
and the factor 1.06 is the value of G(24) . A proper 
derivation of this result will show that a factor 

(ne Lois Dn ) due to the reverse current has also 
been ignored in this intuitive approach to the tunnel 
integral. The major source of discrepancy between the 
observed value of (eee) ahd the above expression 
however will be due to the factor (= V2 a, “This Caeter 
departs seriously from the BCS function G for values of 


NM=A as may be seen from Table A 3.1, 86 that it 


does not represent the distribution B, well in this region, 
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It does not represent the actual Fe very well near the 

gap edge either because of the broadening effects. However, 
the above expression has the advantage that it may be readily 
evaluated since the integral reduces to a gamma function. 


One obtains 


= at (eeC 2x) x a r\/2 -£ 
Tiga) = See = toe (Fr) € 
\ La 2s 
wile face (Ey 2 a : CB aa) 


In this approximation then, one has a simple 
connection between the temperature YT and the observables 
relCZA) and- “A. The observed values of [(2) as a 
function of the observed temperature ora are ~ptovted in 
Figure A 3.2 for two specimens along with equation A 3.1 
as a function of temperature JT . One would expect the 
observed points to at least lie parallel to equation A 3.1, 
since the exponential term due to the energy gap should be 
dominant in any reasonable expression for 1(2M4°) . Some 
experimental results published by Giaever are shown in 
Figure A 3.2 and they confirm this argument. One is forced 
to conclude than that the tunnel junction temperature for 
specimen Al-34 was not the same as the temperature of the 
resistance thermometer with which Tg was observed. Better 
agreement was obtained for specimen Al-42, probably as a 
result of better evacuation of the vacuum chamber of the 


cryostat and a more liberal use of solder at the connection 
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Figure A 3.2 


Graphical presentation of equation A 3.1 along 
with some experimental results. The experimental 
results of Giaever have been obtained from one of 
his publications (Giaever et al, 1962). The 
temperature scale shown has been calculated on 


the basis of Alo)= 0-20 mv. 
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between the He? chamber on which the resistance thermometer 
was mounted and the specimen chamber, but it may be seen 
that further improvement will be necessary if experiments 
are to be conducted which require a knowledge of the tunnel 


junction temperature to within a few millidegrees. 


io “eH sdd meowsed 
sit boas bestavom esw 


yiqmi teadiut gad? 


SS$ZouHnNOS sd os sis 


vyteteqnet aolsonut 


95. 


APPENDIX IV 


DISTORTIONS INTRODUCED BY THE PROBE DISTRIBUTION 


It was stated in Chapter III that the normalized 
dynamic conductance © should be quite similar to the 
normalized density of states in one of the superconductors, 

Peo. , for certain probe distributions P., . The 
correspondence is not complete however, so that © must be 
viewed as a distorted form of jou . The size and shape 
of these distortions will be derived here for some simple 
forms of Po. . 

It has been shown that the variation of the 
tunnelling probablility with applied voltage has a neglibible 
influence on our results, so we will assume a constant 
transition probability here in order to simplify calculations. 
We will also assume that there is a negligible number of 
thermally excited electrons and holes in the distributions 
for the superconductors. The expressions for the tunnel 


current between the superconductors may then be written 


(s(v) = cf Riu) BR. (w-v) dw , 


V 
c [Roo Ryoodx . 
o 


The last expression is obtained by Setting 2s V=t And 
noting that the distributions are symmetric in their 
arguments. The current between these metals when they are 


normal is SG are a . The normalized dynamic 
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conductance is therefore 


Vv 
T= 4 fRovev) Ruwodz - 
; (4.1) 

For cases where My, is normal, we will use 
the Fermi function £(W,T) for ERG in order to take 
into account the effect of a non-zero temperature on the 
probe distribution. 

We now take the point of view generally used by 
electrical engineers and consider equation A 4.1 asa 
mathematical system with one input quantity ne and one 
Output quantity o& - The two quantities are not 
identical due to the distortion introduced by the system 
parameter Ine . We then calculate the distortion 
introduced into some simple forms for Fox . For these 
Simple forms we choose a step function, a ramp function, 
and a parabola, as indicated in Figure A 4.1. The 
voltage origin for R. is chosen at the discontinuity in 
the step function, and the ee functions are taken to be 
neDocor \/ <6 , 80 that higher derivatives of all of them 
are Ciscontinuacus «at V=O . The output results O can 
then be calculated for a given Pe. . The results of some 
calculations of this nature are shown in Figure A 4.1. 
Allowance has been made for the energy gap in choosing the 


origin for the @ results. The energy scale is given 
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Figure A4.1 


Distortions introduced by the probe distributions. 
Parts (a) to (c) are for a superconducting probe 
distribution, while parts (d) and (e) are for a 
normal probe distribution. The distribution Ro 
for parte (a) and (@-is a unit step function; 

for parts (b) and (e) it is a ramp function, and 


for part (c) it @s a parabola. 
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ria te ole (ZN for the superconducting probe distribution 
and in units of YT for the normal probe distribution. The 
unit step function form for R. causes © to be identical 
with the probe distribution, apart from a reversal of 
coordinate direction, so that the shape of the probe 
distribution used for the calculations has not been repeated 
in a separate part of the figure. 

A study of Figure A 4.1 will reveal some of the 
higher derivative results. For example, the derivative 
of the result for the ramp form of Ro is identical with 
the probe distribution, as is the second derivative of the 
parabolic form of Fon . 

The conclusion which we have drawn from the results 
shown in Figure A 4.1 is that the critical point fits shown 
in Chapter V are reasonably valid for energies removed from 


the critical point energy by an amount |X| 2 Ap ‘ 
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